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Abstract. We study the cumulants and their generating functions of the probability 
distributions of the conductance, shot noise and Wigner delay time in ballistic quantum 
\ dots. Our approach is based on the integrable theory of certain matrix integrals and 

applies to all the symmetry classes (3 £ {1,2,4} of Random Matrix Theory. We compute 
the weak localization corrections to the mixed cumulants of the conductance and shot 
noise for /3 = 1,4, thus proving a number of conjectures of Khoruzhenko et al. [49j . We 
derive differential equations that characterize the cumulant generating functions for all 
/3 £ {1,2,4}. Furthermore, we show that the cumulant generating function of the Wigner 
delay time can be expressed in terms of the Painleve III' transcendant. This allows us 
to study properties of the cumulants of the Wigner delay time in the asymptotic limit 
n — > oo. Finally, for all the symmetry classes and for any number of open channels, we 
derive a set of recurrence relations that are very efficient for computing cumulants at all 
' orders. 
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1. Introduction 



1.1. Background. Many successful approaches to tackle problems in quantum transport 
are based on the the Landauer-Biittiker scattering theory of electronic conduction (see, 
for example, the review article |l4J and references therein). It applies to two very general 
types of systems: the first are mesoscopic conductors confined in space, often referred 
to as quantum dots, through which the electric current flows via two point contacts; the 
second are quasi-one dimensional wires containing scattering impurities. The physical 
dimensions of these systems are small enough that the quantum mechanical phase coherence 
of the electron is important and its properties cannot be understood in terms of classical 
mechanics; at the same time they are large enough that a statistical description of the 
electrical current becomes meaningful. 



In 1985 Altshuler [ll| and independently Lee and Stone [56J| discovered that the statistical 
fluctuations of the conductance in disordered wires are universal. More precisely, they are 
independent of the dimensions of the sample or the strength of the disorder provided that 
the aspect ratio length/width ^> 1, the length is much longer than the electron mean free 
path and much shorter than the localization length. Soon afterwards this phenomenon 
was observed experimentally [86]. Approximately at the same time, in a different area 
of physics, it was discovered that as h — > the correlations of the spectra of quantum 
systems with a chaotic classical limit are the same as those of the eigenvalues of random 
matrices from appropriate ensembles 18|, |2l| . It was then realised that Random Matrix 



Theory (RMT) could provide the mathematical framework to develop a statistical theory 
of quantum transport that would account for the universality of the fluctuations of the 
electric current [H, [l3|, 

The main hypothesis behind the RMT approach to quantum transport is that the aver- 
age dwell time of the scattering electron is much larger than the Ehrenfest time, so that 
the system-specific features of the conductor become negligible and physical quantities can 
computed through appropriate ensemble averages, which depend only on the symmetries 
of the system. In this paper we focus on quantum dots connected through the environ- 
ment through ideal leads, or, equivalently, through ballistic point contacts. Under these 
assumptions, the scattering matrix is uniformly distributed in one of Dyson's circular en- 



sembles (12, |13|, ll9|, |20|, |45|: the COE (0 = 1), CUE (p = 2) and CSE {(3 = 4), depending 
on whether the system has time-reversal and/or spin-rotational symmetry. In supercon- 
ducting (Andreev) quantum dots additional constraints imply that the relevant symmetry 
classes are those discovered by Zirnbauer [88| and Altland and Zirnbauer 0, Hoi ]. 
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The subjects of this paper are the cumulants and their generating functions of the 
probability distributions of the electrical conductance, the shot noise, which is the time 
average of the current fluctuations at zero temperature, and the Wigner delay time, which 
is a measure of the extra time an electron spends in the scattering region. The Landauer 
formula expresses the conductance as a sum of the transmission eigenvalues, 7\, . . . ,T n , 
where n is the number of travelling modes (quantum channels) available to the electron 
wavefunction incoming, or outgoing, at one of the two edges of the conductor. It follows 
that both conductance and shot noise are linear statistics of Ti, . . . ,T n . More precisely, 
they are 



(1.1) 



C7/C7 :=E T . and P/Po^Y^TAl-Ti) 



respectively. The constant G := 2e 2 /h is the unit of quantum conductance and P : = 
2eVGo. Here V is the voltage difference between the two edges of the conductor. In what 
follows we shall always choose units where Go = Pq = 1. The Wigner delay time is the 
average of the proper delay times, t±, . . . , r n , namely 



1.2) 



n 

n ^ 



31 



3=1 



where in this case n is the total number of quantum channels. As we shall discuss in detail 
in Sec. I1.2| RMT predicts that the random variables Tj are distributed like the spectra 



of matrices in the Jacobi Ensembles [12] , 1 14 . l28l . l35| . while the proper delay times like the 



inverse of the eigenvalues of matrices in the Laguerre Ensembles [241 ] . In the limit n — > oo 
the conductance, shot noise and Wigner delay time converge to normal random variables 
(see, e.g., 0, Q Q). 

In this article we compute the first corrections to the semiclassical limit n — > oo, often 
referred to as weak localisation corrections, of the mixed cumulants of G and P for — 1,4. 
We also calculate the leading order contributions to the cumulants of the Wigner delay 
time. In addition, we derive a set of differential equations that are satisfied by the cumulant 
generating functions for all G {1,2,4}. In particular, we prove that when = 2 the 
cumulant generating function of rw can be expressed in terms of a solution of the Painleve 
III' equation. 

For over twenty years the universal properties of both linear and non-linear statistics 



of the transmission eigenvalues and t 



using both RMTJ37 J38 



orbit theory [15H17 



23 



49 



41 



54J55U57J59 



5l|-|53|,|58 



re proper d elay times have been extensively studied 
SH, 0, SI BS, H dUl and periodic 



65 



62 



66 



63 



73 



81 



821 ] . Nevertheless, the higher moments 



and cumulants of G, P and r w appeared to be elusive for a long time and considerable 
effort was spent to investigate the lower cumulants 38, 74-7^]. This is hardly surprising, as 
the global fluctuations of the spectra of random matrices are particularly difficult to probe 
and few results are available [22j, |32|, |33j, |47J, |50j. In the case of the Wigner delay time, the 
problem is further complicated by the fact that the right-hand side of Eq. ( 11. 2p is a singular 
linear statistic on the eigenvalues of matrices in the Laguerre Ensembles. Furthermore, the 
appropriate Laguerre Ensemble contains an unusual n-dependent exponent (see Eq. (II. lip ). 
As is often the case in RMT, ensembles with symmetries = 1 and = 4 are more difficult 
to deal with than those with = 2. 



Recently Osipov and Kanzieper [70j, |71| made substantial progress by using the theory of 
integrable systems combined with RMT to study the cumulant generating functions of the 
conductance and shot noise. They proved that when time-reversal symmetry is broken, the 
cumulant generating function of the conductance can be expressed in terms of a solution to 
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the Painleve V equation. They also computed the leading order terms of the cumulants in 
the limit as n — » oo. Independently, Novaes [68| using Selberg's integral discovered a finite 
n formula for the moments of the conductance for (3 = 2 as a sum over partitions. Soon 
afterward Khoruzhenko et al [49] made a further considerable contribution by combining 
generalizations of Selberg's integral with the theory of symmetric functions. In particular, 
these techniques allowed them to compute the first asymptotic correction of the cumulants 
of conductance and shot noise when (3 = 2. They were able to formulate conjectures for 
the same quantities for ensembles with (3 = 1 symmetries, which until then had been much 
more difficult to tackle than the case /3 = 2. In a recent letter Vidal and Kanzieper j83| 
computed the joint probability density function of the reflection eigenvalues in quantum 
dots with broken time reversal and coupled to the environment via point contacts with 
tunnel barriers. 



Our approac 
introduced in 
van Moerbeke 



i is based on the integrable theory of certain classes of matrix integrals first 



40, 



64j . [80j and then developed to its full potential by Adler, Shiota and 
J 4 1 and by Adler and van Moerbeke Using this formalism, we show 

that the mixed cumulant generating function of the conductance and shot noise satisfies 
a certain non-linear PDE, which can be reduced to a set of recurrence relations for the 
cumulants of the conductance and the joint cumulants of conductance and shot noise. By 
looking at the leading order contribution to such difference equations, we can prove the 
conjectures of Khoruzhenko et al. 1491 on the cumulants of conductance and shot noise 
when = 1, as well as perform the same analysis when (3 = 4. 

These techniques are powerful enough to cope with singularities of the linear statis- 
tics (11.21) too. Indeed, we can apply the same approach to the distribution of rw and 
prove a number of finite n properties that were previously beyond reach. These include 
recurrence relations for the cumulants for all (3 G {1,2,4}, as well as the link between the 
Painleve III' equation and the cumulant generating function when (3 = 2. 



1.2. The Landauer-Biittiker Theory and Random Matrix Theory. The dynamics 
of an electron in a mesoscopic conductor can be studied by looking at its scattering by a 
two-dimensional ballistic cavity whose classical dynamics is chaotic. The scattering region 
is connected to two electron reservoirs in equilibrium at zero temperature and Fermi energy 
Ep by two ideal leads, which support a finite number of propagating modes (quantum 
channels). At low voltage the electron-electron interactions are negligible and the scattering 
is elastic. Therefore, the scattering matrix at energy E relative to Ep is unitary and has 
the block structure 

(1.3) S(E) := ( Tmxm V n> ) • 

\t-n xm r nxn / 

Here m and n are the numbers of channels in the left and right leads, respectively; the sub- 
blocks r mxm , t nxm are the reflection and transmission matrices through the left lead and 
r' nxn and t' mxn those through the right lead. Without loss of generality we shall assume 
that m > n. The unitary matrix in (11.31) is often referred to as the Landauer-Biittiker 
scattering matrix. In this setting the semiclassical limit h — > is equivalent to n — > oo 
while the ratio m/n remains finite. 

The transmission eigenvalues T 1; . . . , T n are the eigenvalues of the Hermitian matrix tv 
as well as of t'H'; thus, formulae (11.11) become 



(1.4) 



G = Trtt j! and P = Tr {tt\l - tt ] )) . 
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Since S(E) is unitary, Tj e [0, 1] for j = 1, . . . ,n. The proper delay times are the eigen- 
values of the Wigner-Smith time-delay matrix, which is defined by 

dS(E) 



;i.5) 



Q ■= -ihS-^E) 



dE 



Thus, the Wigner delay time is the normalized trace 
(1.6) 



r w 



-TrQ. 

n 



Note that when we discuss the conductance and shot noise, n denotes the number of 
quantum channels in one lead and not the dimensions of the scattering matrix, which are 
(m + n) x (m + n). When we consider the Wigner delay time, n is the total number of 
channels; thus, in this case the dimensions of S(E) and of Q are n x n. We adopt this 
convention to simplify the notation in the rest of the paper, since n will always be the 
asymptotic parameter. 

Within the RMT approach to quantum transport, the transmission matrix ttf belongs 
to one of the Jacobi Ensembles; therefore, the joint probability density function (j.p.d.f.) 
of Ti, . . . ,T n is 

(1.7) 



12, 14, 28, 35 



Pj(Tx 



T 



1 n 

— TTt c 



:i-t< 



,8/2 



n pi 

l<j<k<n 



Ta 



where a = f (m — n + 1) — 1 and 



'U 



[0.1]' 



J-J^ (1 



T 

± 3 



,5/2 



\Tk — Tj I dTi ■ ■ ■ dT n 



n-l 

n 



a 



n 

+ 1 + 3 (5/2) T {5/2 + 1 + j (3/2) T (l + (j + 1)0/2) 



r(a + 5/2 + 2 + (n + j - 1) /3/2)r (1 + (5/2) 



is Selberg's integral (see, e.g., [36|, Ch. 4). If the scattering matrix S(E) belongs to 
one of Dyson's ensembles then 5 = 0. When additional constraints are imposed new 
symmetry classes may arise. In mesoscopic physics this happens when a chaotic cavity is 
in contact with a superconductor (Andreev quantum dots). The symmetry classes for these 
systems were predicted by Zirnbauer [88[ and Altland and Zirnbauer js), [Iq[. Duenez 31 



studied these symmetry classes too, and extended Zirnbauer's theory. The new ensembles 
modelling the scattering matrix of Andreev quantum dots are the compact symmetric 
spaces that, according to Cartan's classification, are of symmetry types C, CI, D and DHL 
Such symmetries affect the dependence on the integers ((3,5) of the j.p.d.f. (\1.7\i . More 



precisely, the connection between the symmetry class of S(E) and (/3, 5) is the following [28 

((3, 5) = (4, 2), symmetry type C; (/?, 5) = (2, 1), symmetry type CI; 

((3, 5) = (1, —1), symmetry type D; ((3,5) = (2, —1), symmetry type DHL 

Physically, they correspond to different combinations of spin-rotation and time-reversal 
symmetriesEI 

If the Ti,...,T n were independent random variables, then the central limit theorem 
asserts that the variances of G and P would grow in proportion to n; instead, it remains 



1 It is important to emphasise that while /3 = 1 and (3 = 4 label Wigner-Dyson symmetry classes that 
are time-reversal invariant, in the Altland- Zirnbauer classification (9L HoL 31 . 88j j they refer to systems in 
which time-reversal symmetry is broken 1^, 28| . 
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finite in the limit n — > oo. More specifically, we have 12|, [13j, |44j, |45 



1.9) 



lim Var(G) = — 

n->oo 



and lim Var(P) 



64/3 



This phenomenon is a manifestation of the universal conductance fluctuations and is 
due to the strong correlations among the transmission eigenvalues, caused by the fac- 



tor Eli 



<j<k<n 



|Tfc — Tjf in the right-hand side of (11.71) . It is a common property for the 



Similar central limit theorems on linear statistics of eigenval- 

46 . the Gaussian, Jacobi 



29. 3( 



spectra of random matrices. 

ues have been proved for the classical compact groups 
and Laguerre Ensembles [H, 33, 47]. 

Braun et al. [231] and Heusler et al. 41] developed a semiclassical derivation of the average 
and variance of G as well as of the average of P, which at the time was still unknown, 
to all orders in 1/n. Shortly after, Savin and Sommers [76| reproduced the semiclassical 
results using Selberg's integral. Subsequently, the same techniques allowed Sommers et 
al. 79J] and Savin et al. [77] to compute the variance of P and the first four cumulants of 
G non-perturbatively. 

The probability distributions of P and G are strongly non-Gaussian for small numbers 



of quantum channels [491 . |54], |79[ . In the large- n limit they approach a universal Gaussian 



curve and contain weak singularities which have been investigated using large deviation 



estimates by Vivo et al. [84j,|85|]. Furthermore, Khoruzhenko et al. [49] derived exact Fourier 
series representations of these distributions; however, the coefficients of such series involve 
n x n determinants ((3 = 2) or Pfaffians (/3 = 1 or (3 — 4) which are difficult to handle 
explicitly. 

When the scattering matrix belongs to one of Dyson's circular ensembles, then the j.p. d.f. 
of the proper delay times is 24 



;i.io) 



p 



n 



where Z is a normalization constant, 



;i.ii) 



^ n 

j=l ^<j<k<n 

th is the Heisenberg time and 

3(3n 



+ 2-/3. 



1. 



n. 



turns fll.lOp into 



In our choice of units th = n. The substitution Xj = 1/tj, j 
the j.p. d.f. of the eigenvalues of matrices in the Laguerre Ensembles. 

The singularities in the right-hand side of Eq. f ll.lOp make the Wigner delay time con- 
siderably more challenging to study than the conductance and shot noise; indeed, much 
less is known about the probability distribution of Tyj. Lehmann et al. [57| studied the 
parametric correlations of T\y when /3 = 1 for small and large number of quantum channels 
using supersymmetry techniques in RMT. Fyodorov and Sommers 38[ computed the finite 
n parametric correlations as well as the variance when (3 = 2. Furthermore, Fyodorov 



et al. [37j derived the crossover of the same quantities to systems with (3 = 1. Kuipers 
and Sieber 53[ computed such correlations as well as the variance at leading order in n 
using periodic orbit theory for (3 = 1 and (3 = 2 symmetries. Semiclassical calculations to 
all orders in 1/n are still an open problem in this case. The distribution of the Wigner 
time delay for small numbers of channels as well as its transition to non-ideal coupling was 
treated in [73, 25, 78] for each of the Wigner-Dyson symmetry classes. 
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2. Main Results and Discussions 

2.1. Conductance and Shot Noise. The definitions ( II. ip and the j.p.d.f. ( I1.7P imply 

that the random variables G and P are correlated. The jomt moment generating function 

M^(z,w):=^ UT?(1- Tj f> 

(2.1) Cn ' /[0 ' 11 "j=i 

xe3tp(--zT J --«>r j -(l-T i )) JJ m-T/dTx-'-dTn 

l<j<k<n 

contains all the information on the joint probability density function of G and P. By 
definition the joint cumulants nf^ are the coefficients the Taylor expansion 

(2.2) aWfrw) := log^f (*,«,) = EE m 

i=0 fc=0 

Since A^i^(0, 0) = 1, Kqo = 0- The cumulants of the conductance or of the shot noise are 
obtained as special cases by setting w = or z = separately in ( 12.21) ; throughout this 



paper we will reserve the notation for the cumulants of the conductance. 

As n, m — > oo and in — n remains finite, the averages of the conductance and shot noise 
are given by, 

(2.3) = f + and «$ = | + 0(1); 

the limit of the variances as n — > oo are given in Eq. (II .9p . 

One of the main problems addressed in this work concerns the asymptotics of the higher 
cumulants. 

Theorem 2.1. Let = 1,4 and suppose that a is independent of n, i.e. m = n + C for 
some constant C . Then, we have 

(2.4a) lim n l+k K<fi> = (8/2 -a)(a + 5/2 + 2- (3)^ + 



- /3(2/3)*(4/3)* 



\k ■ 

y\^y> K^y. 

for I odd and (I, k) ^ (1, 0), 



(2.4b) lim n'-^ = (-1)* f ^ - l) ( / + * " ^ E ( * + / 1 H HO"', 

1 ; ™ 1 ; V2 ; (4/3) fc (4/?y ^ Vj + w W ; ' 

for I even and (/, k) <£ {(0, 0), (0, 1), (0, 2), (2, 0)}. 

There are some interesting consequences of this theorem that are worth emphasizing. 

Remark 2.2. Let (3 = 1 and 5 = 0. In other words, consider normal (non-superconducting) 
quantum dots with time-reversal and spin-rotation symmetries. Assume also that the 
number of quantum channels is the same in both leads, i.e. m = n and a = —1/2, 
see Eq. ( II ,7p . Khoruzhenko, Savin and Sommers 49] conjectured the weak localization 
corrections of the higher cumulants of the conductance: 

(2.5a) lim n 1 ^ = ^T* , for 1 odd > 



2 l+2 



(2.5b) to = J-L^Q 2 ), for (even 

with I > 2. Equations (I2.5a|) and (I2.5bl) are particular cases of Theorem 12.11 
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In the same article it was also conjectured that at leading order the higher cumulants of 
the shot noise power are 

(2.6a) Jim n k K^ k = for k odd, 



n— >oo 



(2.6b) Ito = -<*Z*>! Q 2 ) , for k even, 

with k > 2. When / = the sum in (12.4bl) is given by 



k 



(2.7) j=0 



2jV*Y 2 )-i = p^-i)(_ 3 )- (- 1 ) fe/ ^ /7T 



j7W v 7 fc v 7 r(-| + |)r(| + i) 




for k odd, 
for /c even, 



where pf: (x) denotes the Jacobi polynomial (see Appendix [C] for the definition). The 
identity (12.71) reduces formula (12.4bj) to the right-hand side of (12.6bj) when k is even; 
however, mutual cancellations occur in the sum (12. 7p when k is odd, which yields 

(2.8) lim n^K^l = 0. 

n— >oo 

This limit implies that the first corrections to the odd cumulants of P are of subleading 
order in 1/n and is consistent with the conjecture (I2.6al) . These cancellations explain the 
unusual "staircase" behaviour in the decay of the cumulants of the shot noise, first observed 



in 



49 



As we shall discuss in Sees. I5.3[ I5.4l and 15. 5[ the last step in proving theorem 12. II consists 
in solving the asymptotic limit of certain difference equations. As it turns out, proving 
formula (12 .6 aft requires computing the next to leading order of the mixed cumulants k^. 
Although in principle our techniques could be carried out at the next order in 1/n, the 
calculations become substantially harder (see Sec. 15. 5p . 

Remark 2.3. Theorem 12 . 1 1 reveals some general qualitative differences among the proba- 
bility distributions of G and P for the various symmetry classes. 

In the special case ((3,5) = (2,0), i.e. when the Landauer-Biittiker scattering matrix 
belongs to the CUE, the cumulants of conductance and shot noise were computed to leading 



order by Osipov and Kanzieper [70|, |7jj. They are of subleading order in 1/n compared to 
formulae (I2.4p , which means that the convergence of G and P to a normal random variable 
is slower when (3 = 1 and (3 = 4. Similarly, the fact that the mixed cumulants are of lower 
order for — 2 j7l[ implies that the correlations between G and P are stronger when (3 = 1 
and (3 = 4. 

For symmetry classes with (3 = 1, the right-hand side of Eq. f!2.4bj) is negative, while it 
is positive when (3 = 4. Thus, in general, in the former case the probability distributions 
have a lower and wider peak around the mean than in the latter. 



Some of the intermediate results leading to Theorem 12.11 are interesting in their own 
right and it is worth discussing the main ideas behind the proof. By suitably deforming 
the integrand of (12.11) with an infinite sequence of "time variables" t = (ti,t 2 , . . .), it can 
be shown that when (3 = 1 or (3 = 4, the resulting integrals are connected to an infinite 
hierarchy of non-linear PDEs known as the Pfaff-KP hierarchy 0]. The first non-trivial 
member of this hierarchy is the Pfaff-KP equation. In addition, the deformations of the 
integrals (12. ip satisfy an infinite set of linear PDEs, known as Virasoro constraints. We 
outline this theory in Sec. 13. 1[ 
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The Pfaff-KP equation and the Virasoro constraints can be used to find a partial differ- 
ential equation that the cumulant generating function (12. 2p satisfies. 



Theorem 2.4. Let (3 e {1,2,4} and 



(2.9) :-- 



ra(n-l) c „-2 c n+2 ■/ o_i 
(„ + l)(„ +2 ) / c (i)^ v P > 



if P = 2, 



(4) (4) 

n c »-i c «+i if 13 = A, 



n+l / (4) 



where is Selberg's integral (11.81) . Furthermore, define the parameters 

[1 i/0e{l,2}, [2 if p = l. 

The mixed cumulant generating function (12.21) satisfies the nonlinear PDE 

(2 11) n W w - — I 6n W iv I — 1 I 2; n 
(2.11) ^to fe4 +b7 7l)4 «;^ fe2 J 

02 (/?) Q W Q2 W) 

= (126f /p>exp (af_ m + ^ _ j^Cffl) . 

We report the explicit formulae for the parameters (I2.9P in Appendix IB. II 

The same approach that leads to (12. lip allows us to derive an ODE for the conductance. 

Theorem 2.5. Denote o~n (z) '■= log Ain^ (z, 0) and let bn , rff\ and i((3) be the same 
coefficients introduced in (12. 9p and (I2.10p . Furthermore, define the polynomials 

(2.12a) p ( i\z) = -~{pn + 2a + 2 - p){pn + 5/2 + a- z), 

(2.12b) pf\z) = (/3n + a- 5/2)z - (f3n + a + 5/2)(/3n + a + 5/2 + 2 - f3), 

(2.12c) pf\z) = -z 3 + 2(/3n + a - 5/2)z 2 

- {pn + a + 5/2 + 6 - 3/3)(/3n + a + 5/2 + 2 - p)z + pz. 

Then 



(2.13) AJS^' + e^ 3 ^") 2 + P{2*>oW» + z{a^f + 4,Vf 

= (12^/pV exp + ^ ^ _ 2t ,Cfl) . 

Remark 2.6. When (/3, 5) = (2, 0) the substitution 

(2.14) a$\z) ^ n(n + a) + za^'(z) 



reduces Eq. (12TT3]) to Painleve V |70|, |Z1| 



Remark 2.7. For P = 1,4, Eqs. (12.111) and (12.131) are differential- difference equations, as 
a n (z) appears in the right-hand side with different values of n. This makes the analysis 
in the cases ft = 1, 4 more difficult. Note that this aspect of the equations disappears when 



10 



F. MEZZADRI AND N. J. SIMM 



fi = 2, as then ft„ = 0. In particular, Eq. (12. lip reduces to the zero-temperature limit of 



a result derived in 71 



Remark 2.8. Since we restrict (i to take only the values 1, 2 and 4, Selberg's integral and 
the coefficient bn^ reduce to a rational function of a, 5 and n (see appendix IB.lj) . 

By introducing the power series expansions of a^P (z, w) and of (z) into the PDE (12.111) 
and the ODE (12. 13p . we obtain two difference equations for the cumulants and nf . 

The K^'s serve as initial conditions of the recurrence relation for the joint cumulants nf^. 
In order to prove theorem 12. II such recurrence relations need to be solved in the asymptotic 
limit n — > oo. This will be achieved in Sees. H] and El 

Why is it necessary to look at the mixed cumulant generating function and not simply 
at those of the conductance and shot noise separately? The main reason is that the shot 
noise is quadratic in the transmission eigenvalues, and consequently, the determination of 
a non-linear ODE satisfied by alf\o,w) is much more complicated than for a^\z,0). In 
a certain special case, we point out that o~n\Q,w) satisfies an ODE related to Painleve V; 
this is discussed in Sec. 12.31 More generally, the derivation of a non-linear ODE satisfied 
by ^(0 , w) is currently an open problem. 

2.2. The Wigner Delay Time. As we discussed in the Introduction, the statistical the- 
ory of the Wigner delay time presents additional difficulties compared with the conduc- 
tance. Nevertheless, the connection with integrable systems is powerful enough for us to 
obtain substantial results. 

Let us set \j = 2tj/ {(in), j = 1, . . . , n, and define 

(2.15) M^(z):=^ m [ TTA- b exp(-^ + ^ TT \X k - \/ dX 1 . . . d\ n , 

where b is the same parameter that appears in (ll.lOp and rrin is a normalization constant 
such that Mn (0) = 1, namely 

r21fi , mW . fr 1 r(l + (j + lWQ - 1 + U - 2n + 2)0/2) 

{2Ab) m " -11 r (i + #2) 

This is a particular limit of Selberg's integral (II ,8p (see, e.g., [Ifl], Sec. 4.7.1). The quantity 
M^t\f3y/2) is the moment generating function of the Wigner delay time (11.101) . 

Theorem 2.9. Let (3 e {1,2,4} and define 

(2.17) £(%):= log Mf) (4 

the parameter 



n(n-l) r4-2 m n+2 
(n+l)(n+2) ( rn m\ 2 



ifP = 2, 

if (3 = 4 



W m (4) 



n+1 (ra^ 



(2.18) dSf) ,- 



and the polynomials 

(2.19a) p { f\z) = 2z-(b-Pn-2 + P)(b-Pn), 

(2.19b) pf \z) = Az 2 + Pz - (b - (in - 6 + 3/3) (b - (in - 2 + p) z. 
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The function ^\z) is a solution of the ODE 

(2.20) V [%V" + 2^"' + 0z (CD 2 + Wz 2 ^"^' 

_ Yldj 3 / m W) _ W) \ 

The quantities vffl and i(/3) were defined in Eq. (I2.10p and b in (II. lip . 

The approach outlined for conductance and shot noise can be applied to the Wigner delay 
time too. However, there are two features of the j.p.d.f. (ll.lOp that increase difficulties of 
the asymptotic analysis substantially. First, the exponent (11. lip is proportional to the 
dimension n, while the parameters a and S that appear in the j.p.d.f. of the transmission 
eigenvalues are taken to be independent of n. Second, the singularity at the origin in the 
integrand of (12.151) implies that Mn(z) is not analytic, therefore only a finite number of 
moments exist. Thus, an infinite power series is replaced by the asymptotic expansion 

(2-21) ^) = £^(|) +o&), z^O, 

where the K^'s are by definition the cumulants of the random variable (11.21) . A brief 
inspection of the integral in the right-hand side of (I2.15P shows that 

q := \b - 2 - /3 (n - 1)J , 

where |_^J denotes the greatest integer such that 

[x\ < x, x G KL 

Because of the increased technical difficulties, our asymptotic analysis focuses principally 
on the simpler case (3 = 2. Our starting point is the following result, an immediate 
consequence of Theorem 12.91 

Corollary 2.10. Let (3 = 2 and set H n (—z) = z^^'(z). Then H n (z) satisfies the non- 
linear second order ODE 

(222) i zH ") 2 = A M(KT - K) 

- (MK) 2 -(4:Z + (b-2n) 2 )H' n -2n(b-2n))H' n + n 2 . 
Proof. Substituting H n (z) into (12.201) gives 

(2.23) z 2 H';' + zH'l + 6z{H' n ) 2 - AH n H' n - (4z + (b - 2n) 2 )H' n + 2H n - n(b - 2n) = 0. 
This equation falls into one of the Chazy classes of ODEs: 

(2-24) L'" + P -L» + ^L>) 2 - ^LL< + P -L 2 + - + ™ = 0, 



where P, Q and R are polynomials. Cosgrove [26|, |27| found a first integral of (12. 24ft : 
(2.25) (L") 2 + ((P(L') 2 + QL' + R)L' - (P'(L') 2 + Q'L' + R')L 

+ \{P"U + Q")L 2 - -P"'l? + C) = 0, 
2 o 
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where C is an integration constant. Setting 

L(z) = H n {z), P{z) = z, Q{z) = -z- [h ~^ n) \ 

R(z) = - n(b - 2n) , C = -n 2 

completes the proof. □ 

Remark 2.11. When ft — 2 Eq. f |2.20j) was studied by Osipov and Kanzieper [69} in the 
context of bosonic replica field theories, who realized that it can be reduced to Painleve 

III . Chen and Its [25] made a detailed study of the partition function ^-MA (—z) and 
showed that 

- 1 - ^ 

is the Jimbo-Miwa-Okamoto a-function for Painleve III'. As far as we are aware, this is 
the first time that the probability distribution of the Wigner delay time is associated to 
Painleve III'. 

When ft = 2 Eq. f)2.20p can be used to obtain a recursion relation for the leading order 
term of the cumulants. Define 

2^-27^(2) 

(2.26) p, = lim 1 , / > 0. 

n->oo [I — L)l 

as well as the generating function 

oo 

(2.27) F(z)=J2p^- 

i=l 

Theorem 2.12. The limit (I2.26j) exists and satisfies the recurrence relation 

2-1 

(2.28) (I + l) Pl+1 = 2(2/ - l) Pl + 2 ^(3z + 1)(Z - i) Pi+1 pi-i, 

i=0 

wrf/i initial condition pi = 1. Furthermore, the generating function F{z) is the following: 

(2.29) F(z) = 3/2 - (3/2)fi(z) + 2zVt(zf + 3^(^) 2 + 2ztt{z), 
where Q(z) admits the power series expansion 

oo 

(2.30) n(z) = l + J2(kZ k , 

k=l 

where 

a4|eV_7 2 )(EC)G 2 -g-)<-3) 



fe-l-i 



Remark 2.13. It is worth pointing out that the coefficients and hence the leading order 
contribution to the cumulants, are integer numbers. It would be interesting to understand 
the physical reasons behind this combinatorial fact. It is not true when ft = 4 (see Table [2]). 
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By equating the first few powers of z in the ODE ( I2.20p . we obtain non-perturbative 
(finite-n) results for the lower order cumulants: 



(2.32a) 
(2.32b) 



K<f> = 


/3n 


Kf = 


(l 


Kf = 





n(2u + fin) 



u\2u + (3){u - 1) (n + l)(n/3-2)' 

4n(2ou + n(3)(u + nP) 
u 3 (2uj + (3)(co - l)(u - 2) 
96 



(2 ' 32c) (n+l)(n + 2)(n/3-2)(n/3-4)' 

where a; = b — 2 — (3(n — 1). Our formula f!2.32bl) for the variance can be extracted 



from earlier results in the quantum transport literature [37J, [38j. When (3 = 2, it has also 
appeared in the context of wireless communications j60[. To the best of our knowledge, 
our formula (I2.32cp for the third cumulant has never appeared in the literature before. 

From Eq. f!2.32bj) it follows immediately that the variance of the Wigner delay time 
satisfies the asymptotic formula 

2, AS) 4 4 I" 1 



(2.33) + ^+0(0 ^oo. 

So far, only the leading term 4//3 in f!2.33j) is known to be consistent with the semiclassical 
theory (53[. We believe it would be important to obtain a semiclassical derivation of f!2.32bj) 
to all orders in 1/n. 

Beyond the third cumulant, the exact results become more complicated. We conclude 
our discussion of the time delay by giving the explicit expressions for the fourth cumulant 
in each symmetry class: 

(2 K m = 96(53n 2 - 68n - 156) 

1 ' 4 " (n-4)(n + l) 2 (n-2) 2 (n + 3)(n + 2)(n-6)' 



(2.34b) K 4 



(-1) 



12(53n 2 + Mn - 39) 



(2.34c) K t 



(n + 2)(n + l) 2 (n - l){2n - 3){2n - l) 2 (n + 3) ' 
12(53n 2 - 77) 



(2) _ 

4 fn2 _ 1 \2(„2 



n 2_ l)2( n 2_4)( n 2_9)' 



2.3. Shot Noise when n = m and /3 = 2. We end the overview of our results by 
discussing two interesting properties of the shot noise for symmetry classes with — 2 
that have remained unnoticed in the quantum transport literature. These are immediate 
consequences of the definition of the moment generating function fl2.ll) and of previous 



results from Witte et al. 87 and Forrester 34 



When the two scattering leads in the cavity support the same number of quantum 
channels, i.e. n = m, the description of shot noise simplifies considerably if (3 = 2 because 
a = 0. We first discuss how the moment generating function can be expressed in terms 
of the Painleve V transcendent; then, we point out an interesting identity between the 
probability distributions of the shot noise and the conductance. 

2.3.1. Shot Noise and Painleve V. It is well known that gap probabilities in ensembles 
with (3 = 2 are related to Painleve transcendents. In particular for the Gaussian Unitary 
Ensemble we have the following. 
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Proposition 2.14 (Witte, Forrester and Cosgrove [87]). Consider the Gaussian Unitary 
Ensemble of random matrices and let E 2 (w) be the probability that the interval (— oo, —w)U 
(w, oo) contains no eigenvalues. Then, 



(2.35) 

satisfies the non-linear ODE 

(2.36) wR" + 2R! = 2w(w - h) - 2h v / (R + wR') 2 - Aw 2 (w - h)R - 2nw 2 (s - h) 2 , 
where 

(2.37) h = Vw 2 - 2R'. 

The function R(w) can be expressed in terms of Painleve V transcendents with appro- 
priate boundary conditions [871 ]. 

It turns out that the gap probability Esi(w) is directly related to the moment generating 
function of the shot noise. More precisely, we have 



(2.38) 



E 2 {w) = Ce- w2n (2w) n2 M£\0, -Aw 2 ), 



where C is a multiplicative constant that does not affect the definition f l2.35p . The iden- 
tity dZ3g) follows from 



M 2) (o, 



-w) 



(2.39) 



n J[0,l\ n j= i l<j<k<n 



C. 

e wn/4: w -n 2 /2 
(2) 



_rp2 

e j 



) \ \T k -Tj\ 2 dTx-'-dTn 

l<j<k<n 



which is obtained by completing the square in the exponential and by changing variables. 



2.3.2. Shot Noise and Conductance. There is fascinating relation between the distribu- 
tions of the shot noise and the conductance in non-superconductive and superconductive 
quantum dots, respectively. 

Proposition 2.15. Let Pn 2 '°\ and G^ denote the random variables defined 



in (jl.4j) . where 
(2.40) 



JV i: = 



71+ 1 



and No :- 



The subscripts refer to the number of variables in the j.p.d.f. (11. 7p . while the superscripts 
emphasize its dependence on the parameters (j3, 5). In addition we assume that a = 0. We 
have 



(2.41) 



r (2,-l) , ^(2,1) 
p(2,0) A. T ^N 2 



where the notation = denotes equivalence in distribution. 



Proof. The proof follows from a result of Forrester [34J , who exploited the symmetry of 
the integration interval in the right-hand side of (I2.39p . together with the evenness of the 
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Gaussian weight. He thus obtained the identity 

AT- 



M^(0,-w) = e --- TT 27^-^/4 jj \ Tk - Tj \ 2 dT x ...dT Nl 



(2 42) ° n °' 1]Nl 3=1 1 ' ' X 



, ^2 



[0,1)^2 j . =1 l<j<k<N 2 

Changing the variables 1} — > 1 — Tj, j = 1, . . . , n, and w — > 4w gives 



(2.43) MS ) (0,-4w)=MP 1 (- 



-w 



8=-l 



8=1 



Finally, Eq. (I2.4ip follows from the convolution theorem. □ 

When (/3, 5) = (2, —1), then the compact symmetric space of the scattering matrix S is 
of symmetry type Dili; this means that S is orthogonal and self-dual. When (/3, S) = (2, 1) 
the symmetric space is of symmetry type CI; in this case 5* is symplectic and symmetric. 
Note that matrices of symmetry type Dili are real, while matrices of symmetry type CI 
have quaternion elements. In both instances time-reversal is preserved. The Landauer- 
Buttiker scattering matrix S of the system on the right-hand side belongs to the CUE; 
thus, interestingly, time-reversal is broken on the left-hand side of Eq. (I2.4ip but not on 
the right-hand side. 

From a computational perspective, we note that the right-hand side of (12.411) is far more 
amenable to calculation than the left-hand side. As we shall see in Sec. HI the cumulants 
of G^' ^ and G^ can be determined from the one-dimensional difference equation (14.61) . 
while to obtain the cumulants of the shot noise we have to study the partial recurrence 
relation ffl~T6]) . 

3. Integrable Systems and Moment Generating Functions 

We can associate to a scattering observable, such as the Landauer conductance or the 
Wigner delay time, a moment generating function whose integral representation is a parti- 
tion function of one of the orthogonal, unitary or symplectic ensembles of random matrices, 
depending on whether (3 = 1, /3 = 2 or /3 = 4. In this section we describe the general the- 
ory that provides a link between such partition functions and certain integrable non-linear 
differential equations. 

3.1. Matrix Integrals, r-Functions and Integrable Hierarchies. Fundamental to 
our approach is the connection between matrix integrals and exactly solvable models, as 
established in a variety of papers, e.g., (2], 0, 0, H, 64, [80[. These works focus principally 



on the case /3 = 2 and obtain relations to the Toda lattice and KP hierarchies. In contrast, 
the existence of integrable structures related to the (3 = 1,4 matrix integrals 0, @, S] has 
received comparatively little attention. 

The starting point for establishing the relationship with integrable hierarchies is the 
deformed integral 

(3.1) (t) := — I TTp(xj)exp ( V^i4 ) FT \x k - Xjf dx x ■ ■ ■ dx n , 

nl J[A,B]- ]=1 W J l<j<k<n 

defined by introducing an infinite sequence of time variables t = (tx, £2, £3, • • •) into the inte- 
grands of Eqs. (12.11) and (I2.15p . For the mixed moment generating function of conductance 
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and shot noise, [A, B] = [0, 1] and the weight function is 

(3.2) Pz,w(T) = T a (l - T ) s / 2 e~ zT - wT ^- T \ 

For the Wigner delay time, we have [A, B] = [0, oo) and 



(3.3) p z (A) = A- ft exp^-- + ^Aj , 

where z < 0. In the theory of integrable systems the integral ( 13. ip is known as a r-function. 

When (3 = 2 the integral (13 .ip simplifies considerably and can be written as a Hankel 
determinant, which can be studied using orthogonal polynomials. We have 

(3.4) r^(t) = det ( [ B x i+j p(x)e^ Uxl dx) 

\JA J 0<i,j<n-l 

When (3 = 1, however, no such nice structure exists. In this case, for even n the r- 
function (13. ip becomes the Pfaffian of a skew symmetric matrix (the square root of its 
determinant). Namely, (see, e.g., |61j ) 

(3.5) r£\t) = ¥{(\ [ [ y k x l sgn(y-x)e^ u(yt+xl) p(y)p(x)dydx _ 

' 0<k,l<n-l 

Although we shall restrict our attention to the weights defined in ( 13 .2p and (I3.3p . the 
identity ( 13. 5 p holds independently of these choices or the deformation parameters t. For- 
mula ( 13. 5 p is based on an identity of de Bruijn for integrating determinants: 



(3.6) — / det(F i (y 1 ),Gi(y l ),...,F i (y n ),G i (y n )) 0<i<2n _ 1 dyi---dy n 

'[A,B] n 

-B 



111 



Pf (7 {Gi{y)Fj{y) — Fi{y)Gj{y))dy 



0<ij'<2n-l 



where F{ and Gj, i = 0, . . . , n — 1, are two sets of integrable functions. In the left-hand 
side of (I3.6p . the commas separate columns and the index i labels the rows. 

By studying the evolution of the moment matrix (13.51) in the extended t-space, Adler 
and van Moerbeke 0, @, II] discovered a bi-linear identity that applies to r-functions that 
can be expressed as Pfaffians: 

(3.7) / r (i)(t-[ r i]) T « 2 (t'+[r 1 ])e^^-^r-^ a ^: 

= - £ r«(f - [e])r« 2 (t + [Z))eZ"M-m-^- m ^ 

where n and m are even, and or Co are small circles in the complex plane enclosing 
the point £ = oo and £ = 0, respectively. The notation [£] refers to the infinite vector 
(£,£ 2 /2,£ 3 /3,...). 

The identity (13.71) can be used to generate an infinite sequence of integrable hierarchies 
of PDEs involving r„ (t). Evaluating the residues on both sides of (13. 7p and expanding 
up to terms linear in (t — t'), one obtains a hierarchy of PDEs known as the Pfaff-KP 
hierarchy; its first non-trivial member is the Pfaff-KP equation: 

/ r) A r) 2 r) 2 \ / r) 2 \ 2 r (1) r (1) 



dtf dt 2 dtidt 3 J 6 n \dt 2 6 n J ( r W) 2 
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When Eq. (13. 8p was introduced in RMT, it was used to characterize gap formation prob- 
abilities in orthogonal and symplectic ensembles |6|. To the best of our knowledge, this 
article is the first to apply ( 13. 8 p in a physical context for over a decade. 

A similar discussion holds when (3 = 4. One can show that ri, (2t) satisfies the same 
equation as Tn\t), except for one small difference: in the right-hand side the sub-indices 
n ± 2 are replaced by n ± 1 and there is no parity restriction on n. 

Although our emphasis is on the — 1,4 ensembles, we point out that when (3 = 2 
there is a simpler and more widely known analogue of (I3.8P usually referred to as the KP 
equation: 

< 3 - 9 > (| + 3 |- 4 l^) logT » 2) + 6 (| lo ^» 2, ) 2 -°- 

which is the first member of the KP hierarchy. 

It is worth emphasizing that the right-hand side of Eq. (I3.8P is not zero as well as 
containing (t) (/? = 1) or T^ 1 (2t) {(3 = 4). This constitutes an additional complication 
compared to the symmetry classes with (3 = 2. 

Equation (13. 8p is the starting point to obtain the differential equations in Theorems I2.4[ 
12.51 and 12.91 in turn such differential equations lead to recurrence relations for the cumu- 
lants. Deriving these differential equations requires a combination of techniques pioneered 
by Adler and van Moerbeke 0]. More precisely we need extra information on the proper- 
ties of r-functions that is provided by the fact that (13.11) satisfies the Virasoro constraints. 
These are an infinite sequence of linear differential equations with a specific algebraic 
structure, which we discuss next. 



3.2. Virasoro Constraints and /3-Integrals. The Virasoro constraints were first ap- 
plied to unitary matrix models ((3 = 2) in the context of string theory jio|, 64[. In the 



unitary case, Osipov and Kanzieper [71J were the first to apply these techniques to quan- 
turn transport and used (13.91) to derive a PDE satisfied by log M. n (z, w). The fact that the 
Virasoro constraints apply for general (3 > appears not to have been applied to problems 
in quantum transport until now. 

The integrals (13. ip are well-defined for any real (3 > 0. Usually, when we want to 
emphasize this property, we refer to them as (3 -integrals. The Virasoro constraints we need 
follow from a particular change of integration variables. More generally, let f(x) and g(x 
be real analytic in a neighbourhood of the origin and supported on an interval [A, B] C ] 
In addition impose the boundary conditions 

(3.10) V\x) :=^= ffi M ! , lim f(x)p(x)x k = 0, for all k > 0. 

/ i X ) Li=0 a i X X ~> A > B 

The r-functions (13.11) belong to the class of /3-integrals with weights of the form 

(3.11) p(x) := e~ vix) , xe[A,B]CR. 

It is straightforward to check that the boundary conditions ( I3.10P are satisfied by both 
weights (Q and ( 133|) . 



In general the domain of integration of (|3.1I) is E n , where E is a disjoint union of intervals contained 
in [A, B\. In the cases we study, E = [A, B], which slightly simplifies the analysis. 
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Let us introduce the differential operators 

k > 0, 



J, 



(i) 



d 
dt k 



k = 0. 



fc-i 



5 2 



i=l 



2 ™ g 



where t%, 
(3.12a) 

(3.12b) 



, tk are the first k elements of the infinite vector t = (^1,^2,^3 



Define 



+ U + (t + l)[l-^))jf>+n|(»-l)S + l)*M 



(3 



as well as the Virasoro operators 

00 

(3-13) V k := £ (a,4 2 i, n - h 



,t(1) 



k > -1, 



i=0 



where the parameters and 6j, i = 0, 1, 2, ... , are the coefficients of the Taylor expan- 
sion of f(x) and g(x) defined in Eq. (I3.10p . Imposing the invariance of (13. ip under the 
transformation 



(3.14) xj -> xj + ef(xj)xy \ j = 1, . . . , n. 
one can show that (see Appendix 5 of 0]) 

(3.15) V fc T„(t) = 0, k > -1. 

We shall refer to these identities as Virasoro constraints. A more abstract proof of 
Eq. (I3.15P is based on the fact that the /3-integrals (13. ip are the fixed points of vertex 
operators jfj|. 

It is worth emphasizing that the differential operators (I3.12p form an algebra, although 

this fact will not enter directly in our discussion. For all n G Z + the operators S^' n and 
(2) 

J k n obey a Virasoro and Heisenberg algebra 

k 



j(i) j(D 



3 



:5, 



T (2) T (l)' 



j(2) ^(2) 



fc,— 6 



(k £) 



, k 3 -k\ . 

'• I — ^ ) <h: 



with central charge 







1/2 



-1/2^ 



and c' = — . 

(3 2 



r-FUNCTION THEORY OF QUANTUM CHAOTIC TRANSPORT WITH /3 = 1,2,4 



19 



3.3. Virasoro Constraints for Joint Conductance and Shot Noise. To put the 

theory of the matrix integrals described in Sees. 13.11 and 13.21 into our context, let us look 
at the r- function obtained by deforming the joint moment generating function (12.11) . 

Since only linear and quadratic terms in the transmission eigenvalues Tj appear in the 
weight (13.21) . we only need the first two Virasoro constraints. The logarithmic derivative 
of the weight (J32D 

f(T) d 

_a-(a + 5/2 + z + w)T + (z + 3w)T 2 - 2wT 3 
~ T(T - 1) 

gives the coefficients of the Taylor expansions of f(T) and g{T): 

ao = 0, ai = —1, a 2 = 1, «3 = 0, 

60 — a, bi = —(a + 5/2 + z + w), bi = z + 3w, 63 = — 2w, 
^ = 0, bi = if i > 3. 

The first two Virasoro operators are: 

V- = £ «< - S-) + (° < + * + ^2 + w + M + 2 - fl)^ 

2=1 

- (z + 3w) — + 2u> — - fan + n((n - l)/3/2 + 1) 



(3.16a 



and 



^=E*(a^-^r)4S-(^ +2 -« + -)lr 



(3.16b) 



+ (a + 2 + 5/2 + w + (n/3 + 3 - 3/3/2)) — 



, 9 d 



The Virasoro operators associated to the deformation of the moment generating function 
of the conductance are obtained simply by setting w = in Eqs. f )3.16al) and fl3.16bj) . 



4. Conductance and Shot Noise for Finite-™ 

The Virasoro constraints allow us to evaluate the Pfaff-KP equation (13. 8p and the KP 
equation (13. 9p at t = 0. This gives the differential equations ( 12. lip and (12 . 13 j) . thus 
proving Theorems 12.41 and 12.51 Then, by inserting the series expansions (12. 2p into these 
differential equations we obtain recurrence relations, which allow both non-perturbative 
and asymptotic analysis of the cumulants as n — > oo. 

4.1. Differential-Difference Equations: the Conductance. As discussed in Sec 12.11 
the cumulants of the conductance provide the initial conditions for the recurrence rela- 
tion for the mixed cumulants k\^. Therefore, we shall start with the proof of Theorem 12.51 
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Let us write down the r-function associated to the weight (13. 2p explicitly 



(4.1) 



i r 11 

r n (t;z,w):=- / H^l-Trf" 
nl Jm n j=1 



x (>xp j -zTj - wT d (l - Tj) +^T}ti \[ \T k - TjfdTi ■ ■ ■ T n . 

i=l J l<j<k<n 



We now study the r-function r n (t; z) := r n (t; z, 0) as well as the moment generating func- 
tion 



(4.2) 



Cn 



t=0 



First note that the definition (14.11) and (13. ip imply the relation 

(4.3) ^Wfrz) = -^frz). 



Eqs. ( 14 .2p and (14. 3p provide the projections at t = of the right-hand side of ( 13. 8 p (for 
both (3 = 1 and = 4) as well as of the partial derivatives with respect to ti in the left- 
hand side. Therefore, in order to complete the proof of Theorem 12.51 we need to express 
the partial derivatives 



(4.4) 



and ^ 

t=0 <%i<9£ 3 



logrf (t; 



t=o 



in terms of the derivatives of ai 0) (z) = log M^\z). This is a standard calculation involving 
the Virasoro operators (13. 16a|) and (13. 16b|) . whose details are provided in Appendix 1X1 



4.2. Non-Perturbative Recurrence Relation for the Cumulants of the Conduc- 
tance. The ODE (I2.13P combined with the Taylor expansion 



(4.5) 



1 z l 



1=1 



provides a non-perturbative recurrence relation that the cumulants satisfy. In Sec. |5] 
we shall study its asymptotic limit as n — > oo. 

Lemma 4.1. The cumulants satisfy the recurrence relation 



i-i 



i=0 



(4.6) A(1)k^ + nfl ( i) 4%K?Mh I) - 1(21 - l)(a - 5/2 + /3n)*f 

- i(i 2)4% = {mwim -w- 2)^1 i > 3 
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(4.7a) 
(4.7b) 



.09) 



n 



(a + 1 + g(w - l)/2) 
a + 5/2 + 2 + /3(n- 1)' 

1 n(2a + 2 + /3(n- 1)) (5 + 2 + /3(n - 1)) 

4 (« + 5/2 + 2 + /3(n - l)) 2 (a + 5/2 + 3 + /3(n - 1)) 



x 



(4.7c) 



K 



(5 + 2a + 4 + (3{n-2)) 
(2a + 5 + 4 + /3(2n-3))' 

2«5 8) ((J/2-a + 2/3/c? ) - /9n) 



(a + 5/2 + 4 + /3(n - 1)) (a + 5/2 + 2 + /?(n - 2)) ' 

The parameter bn^ was defined in Eq. (12.91) and the coefficients A(l) and B(i,l) are 
polynomials in i and I: 

(4.8a) A{1) := r$l{l - - 2) + (31(21 - 1) - (6 - 3(3){a + 5/2 + 0n + 2 - (3)1 
- (a + 5/2 + (3n)(a + 5/2 + (3n + 2 -/?)(/ + 1), 

(4.8b) /) := (2 - X [%1 - i) ((/ - i)(6i - 2) + 3) + ( X g - 1)(J - z) 2 (6* + 2), 

where 

'\ if (3 

2 if/3 

The coefficients /^ , I = 0, 1, 2 . . . are the "reduced moments" , which can be obtained from 
the recurrence relation 



A 1.2 



1,4, 



2. 



(4.9) 

where 
(4.10) 



(0) 



£ 

i=o 



/ - 1 

j 



l = 1,2,... 



,(£)._ , W) 



( 



"2 



ra— >n— i(/3) 



+ AC 



(0) 



n— >n+i{0) 



2k 



and by convention we have set [Iq — 1. TTie integer i((3) was introduced in Eq. (I2.10p . 

This lemma can be proved by inserting the expansion H4.5|) into the ODE (I2.13P and 
performing straightforward, although tedious, calculations. The initial conditions (14. 7p are 
obtained by equating the coefficients of the first three powers of z on both sides of the 
ODE. The right-hand side of Eq. (12 . 13[) contributes only for powers z l such that I > 3 and 
is expanded using the identity 



(4.11) 



exp 



00 03) / 
Pi* 



J=l 



1=1 



l\ 



Since the coefficient /^[_ 3 on the right-hand side of (14. 6 j) involves only the first / — 3 
reduced cumulants, the difference equation (14. 6 j) and the initial conditions (14. 7p define the 
cumulants of conductance for systems with (3 = 1,4 uniquely. It also presents an efficient 
algorithm for their computation. 



Remark 4.2. Formulae (14 .7p generalise previous results |49|, |76|, |77J, |79j to the symmetry 
classes 5^0. The right-hand side of the differential equation (I2.13P does not contribute to 
the recurrence relation until one equates the coefficients of the third power of z. In other 
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words, it does not contribute to the first three cumulants. It is possible to verify using only 
the Virasoro constraints that the formulae in (14.71) are valid for arbitrary > 0. 

Remark 4.3. By setting I = 3 in (14. 6 p we obtain exact formulae for the fourth cumulant 
of conductance in terms of the first three. For example, when /3 = 1, 2 we have 



(4.12a) k 



(i) 

4 



(4.12b) k 



(2) 



-2k { 2 ] + Wk^k^ + 22(4 1) ) 2 + 5^(5/2 - a - n) - 24bi 
(a + 5/2 + n + 4) (4a + 25 + An - 3) 

-2/4 2) + 20/4 2) 4 2) + 32(4 2) ) 2 + 54 2) (<V 2 -a-2n) 
(a + 5/2 + 2n) 2 - 9 



When /3 = 1, formula (14. 12 aft involves the constant bn , whose value is given explicitly 
in Eq. (IB.laj) . At first sight it seems that both cumulants fl4.12aj) and (14.12bj) should be 
of order 0{n~ 2 ) as n — > oo, since Kj is 0(n 2_J ) for j = 1,2 and k% is 0(n~ 3 ). A 
closer inspection reveals that (14.12bj) contains mutually cancelling terms which lead to k± 
being 0(n~ 4 ). When (3 = 1 these cancellations are suppressed by the constant b$ causing 
= 0(n~ 3 ). From the recursion relation (14. 6p . one easily sees that this "staircase" effect 
persists in all higher cumulants of the conductance. When (3 = 4 we observe a similar 
behaviour to the case (3 = 1. 

4.3. Differential-Difference Equation: Joint Conductance and Shot Noise. The 

aim of the this subsection is to prove Theorem 12.41 by deriving the non-linear PDE (12.111) . 
The strategy, as in the proof of Theorem 12.51 is to combine the Pfaff-KP equation (I3.8P 
satisfied by the r-function (14.11) with the Virasoro constraints. 
First observe the trivial identities 

M 

(4.13a) r^(t;z,w) = ^-M^(z, w), 

t=o nl 

(4.13b) ^\t;z,w) = - J^rf )(t;z, W ), 

(4.13c) 8- T W(f,z )W ) = ^-rjf\t;z,w) - ^ (t; z,w). 

Therefore, the partial derivatives of any order with respect to t\ and t% at arbitrary order 
of Tn (t; z, w) can be expressed in terms of the partial derivatives of AiT\z, w). It follows 
that the only missing term in the Pfaff-KP equation (13. 8 p is 

(4-14) — -— logx^ (t,^,w;) 



dhdU 



t=o 



Writing the constraint (13. 16a|) in terms of logr n (t, z, w) and taking the derivative with 
respect to ti, the projection at t = yields the relation 



|__|_ + (a+V2+z+ro+ „ /J+2 _ /J) g 



Q2 



= 0. 

t=o 



Using the relations ( I4.13P we can solve this equation for ( I4.14p . which can then be written 
in terms of t 
PDE HZIIJ) . 



in terms of the derivatives of Ain(z,w). Eventually, after some algebra, we obtain the 
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4.4. Non-Perturbative Recurrence Relation for the Joint Cumulants. The last 
result that we need in order to pursue an asymptotic analysis of the joint cumulants of 
conductance and shot noise is the following. 

Lemma 4.4. The cumulants of the joint probability density function of conductance and 
shot noise obey the recurrence relation 

(4.16) M^SUi - 4%, k -i) ~ 2(a + 5/2 + n/3 + 2 - P)k\% + (21 + 3k + 2)k^ +1 

fc - x /U i\ 1 /A ... ... n#) 



+ Dr ?l,4 /C 2^\ v j ^ ^ / S'+2,i K /-j+2,fc-i-l - ^ 



i=0 x ' j=0 

oo / ,„i\O0 



with boundary conditions ^K^fj = y^fj > w hich are the solutions of the difference 
equation (14. 6p . The coefficients //[^ are defined by the recurrence relation 



<«■") "U=fi(*7 1 )?G) ra 



—j,k—ir u j,i 



i=0 v 7 i=0 

OO / /„ s \ OO 



loiZ/j boundary conditions ^/-tfo) = (^z^) 5wen ^ e solution of Eq. (j4.9j) . TTie 
coefficients rj^ are the "reduced joint cumulants" 



(4.18) r£W« 



n— »n— j(/3) 



T7ie parameters rff\, i(p) and bffl are the same that appeared in Theorem 12.41 and were 
defined in Eqs. ([Z2]) and fl2TTD|) . 

The proof is analogous to that of Lemma 14.11 provided the expansion (14.51) is replaced 
by ([22D, the ODE f[2T3l) by the PDE fl2TTT|) and Eq. flCTl) by 

(OO OO W I A OO OO (P) I k 

EEiir = 1+ EEiir- 
1=0 k=0 / 1=1 k=0 

Remark 4.5. Setting k = in (I4.16P leads to the relation 
(4.19) (I + l)^? = (a + 5/2 + nP + 2 - P)k$ 2>0 , 

which holds for G {1,2,4}. For example, setting I = and using (I4.7bj) gives an exact 
result for the average shot noise: 

m _ 1 n(2a + 2 + P(n - 1)) (5 + 2 + P(n - 1)) (2a + 5 + 4 + g(n - 2)) 
^ ~ 4 (a + 5/2 + 2 + Pin - 1)) (a + 5/2 + 3 + - 1)) (2a + 5 + 4 + /3(2n - 3)) ' 

Instead, setting / = 1 shows that the covariance between conductance and shot noise is 
related in a very simple way to the third cumulant of the conductance given in Eq . (I4.7cp . 



When P = 2 and a = 5 = formula (14. 19[) was derived for all Z > 1 in [71| . Since 
when k = the right-hand side of the recurrence relation (14.161) does not contribute, we 
conjecture that (14.191) holds not only for P G {1, 2, 4}, but also for any real P > 0. 

Remark 4.6. For convenience, in Table [T] we list the first few values of the sequence 
^/iffc) appearing on the right hand side of (I4.6p . They are expressed in terms of the 

reduced joint cumulants r^, defined in Eq. (I4.18p . For ease of notation we have dropped 
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the superscript from the formulae in the table. Interchanging I and k in /i^ leads to 
the same expressions but with the indices of interchanged. 



TABLE 1. The first few mixed reduced cumulants 



1 k 


Hl,k 




1 

1 1 

2 
2 1 


1 

r i,o 
fo,i r i,o + ri,i 
r 2,o + rf fi 
7*2,1 + 2ri j0 r lil + r 0A r 2fi + r ,ir 2 



Remark 4.7. By setting (I, k) = (0, 1) in (I4.16P we obtain a non-perturbative formula for 
the variance of shot noise. When 5 = this quantity was already calculated in 49|, |55|, [77 . 



Our formulae below are also valid for 5^0, and are expressed in terms of formulae (14 .7p 
and (1412]) : 



K (4) _ I 

K °< 2 " 5 



(2/3)(a + 5/2 + An - 2) 2 nf - 4^ 4) - 24(4 4) ) 2 + «<*> + 36«) 
i ((2/3)(a + 5/2 + n + l) 2 /^ - «4 J - + # + 126^ 



5. Asymptotic Analysis 



We shall now prove Theorem 12.11 This means solving the recurrence relations (14. 6|) 
and (14.161) at leading order as n — > oo. It turns out that in this limit, and for each 
(3 G {1,2,4}, Eqs. (14.61) and (14.161) become linear homogeneous recurrence relations which 
can be solved with elementary methods. 

5.1. Cumulants of the Conductance. Firstly, we need to study the asymptotic limit as 
n — > oo of the difference equation (14. 6p . as its solution will provide the boundary conditions 
for the leading order contribution to Eq. (14.161) . 

Lemma 5.1. Let (5 = 1,4 and suppose that a is independent of n, i.e. m = n + C for some 
constant C , where m and n are the quantum channels in the left and right leads. Then, 
the cumulants of the conductance admit the asymptotic expansion 

oo 

(5.1) rep 3 ~ n~ u( - l) Kf\q)n- q , n -)• oo, I > 3, 

where 

if I is even, 
if I is odd. 



(5.2) u{l) = I - e(l) with e{l) 



Furthermore, the leading order contribution k\ (0) satisfies the difference equation 
(5.3) - /3 2 (Z + 1)^(0) + - 1)(4Z - 5)^(0) 

3 • ~ >^ • >^ m 



16/3 ; 



:/(/- l)(/-2)(/- 3 + e(/))(/ -4 + e(/))4! j 3 (0), I > 6, 
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with initial conditions 

(5.4a) 4"'(0) = W2-«)(° + f/2 + 2-g) 

H(A/9 — rv\(, 

(5.4b) «W(0) ( , 

for the odd cumulants and 



(£W _ 3 ( 5 /2 - a) (a + 5/2 + 2 - /3) 



/or i/je even cumulants. 

Since the parameters a and 5 are independent of n, it follows from the finite- n recurrence 
relation ( 14. 6 p and from the formulae for its coefficients A(l), B(i, I) and bn \ see Eqs. (14. 8 p 
and (IB.ip . that the cumulants are rational functions of n; therefore, they admit an asymp- 
totic expansion of the form (15. ip . The initial conditions ( 15. 4 p and (I5.5P are derived using 
the finite- n recurrence relation (14. 6p and by computing the leading order term as n — > oo. 
They are consistent with the expansion (15.11) . 

We shall prove the leading order difference equation (15.31) by appropriately rescaling the 
recurrence (14. 6 p and using a dominant balance argument. The explicit expression (15. 2p for 
the exponent v(l) will follow by induction on Z; that is, we will assume that df\ • • • , d-P 

have the asymptotic expansion (15.11) and prove that does too for any I > 6. We shall 
discuss the left- and right-hand sides of Eq. ( I4.6P separately. 

5.1.1. Proof of the left-hand side of Eq. ( 15. 3p . Recall that 

A(l) := V [^l(l - l)(l - 2) + 01(21 - 1) - (6 - 3/3) (a + 5/2 + 0n + 2 - 0)1 
- (a + 5/2 + 0n)(a + 5/2 + 0n + 2 - 0)(l + 1) 

and 

B(i, I) := (2 - X S3)(J - i) ((/ - i)(6i - 2) + 3) + ( X g - - *) 2 (6* + 2)- 
The second term in the left-hand side of Eq. (15. 3p is 

(5.6) lim n l+ ^~ 2 ( V fl V ( l ) A^B(i, I) - 1(21 - l)(a - 5/2 + 0n)n^ 




- 1(1 -l)(l- 2)<; ) = \l(l - 1)(4Z - 5)^(0), / > 6. 

The computation of this limit is based on the following considerations. Because of the 
asymptotic expansion (15. ip . we see that each term of the sum in Eq. (15. 6p is of order 
^i+i^i— i = 0(n~ l ~ e ^ +1 ) and so it does not contribute at leading order. The exceptions to 
this come only from terms with i G {0, 1, / — 2, / — 1}, which yield systematic contributions 
to the right-hand side of (15. 6p . By taking only these values of % in (15. 6p . together with the 
asymptotic formulae, 

(5.7a) «<0 = 2 + 2^ + O ( n -i), 

(5.7b) K f = ^ + 0(n- 1 ), n^oo, 
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one easily derives ( 15. 6p after a number of delicate cancellations. 

5.1.2. Proof of the right-hand side of Eq. (15. 3p . The right-hand side of the difference equa- 
tion (15.31) is inherited from the right-hand side of the Pfaff-KP equation (13.81) . which is the 
cause of the increased complexity of the symmetry classes with /3 = 1, 4 compared to those 
with (3 = 2. It is not surprising that it requires more attention than the left-hand side. 

The main problem when trying to solve directly the ODE (12.131) or the finite-n recur- 
rence (14.61) is that both are difference equations in n too. However, one would expect that 
this should not affect the leading order asymptotics of Eq. (14.61) as n — )■ oo. Thus, the first 
step is to compute the leading order term of the reduced moments , which by definition 
include cumulants of dimensions n ± i((3). 

Recall that the ^ s are expressed in terms of the reduced cumulants 

(5.8) rf := k (/3) 



via the recurrence relation (14.91) . Since the parameters a and 5 are independent of n, by 
assuming (15.11) and inserting it into (15.81) . we see that admits the asymptotic expansion 

oo 

(5.9) r\ P) ~ n-^ l) - 2 ^ r\ P) (q)n- q , n -> oo, I > 6 

g=0 

and that the leading order term is 

= z(/3) 2 (/- e (/))(/ + l- e (/))^ ) (0). 

In order to relate the asymptotic limit of the reduced moments //[ to that of the r| s, 
it is convenient to use the formula 

rather than the recurrence relation (14.91) . In the right-hand side of (15. lip , tt = {B\, . . . , B p }, 
1 < p < I runs through all the partitions of a set U = {1,2,...,/} of / elements into 
disjoint subsets. The notation \Bi\ indicates the cardinality of the subset Bi C U . Now, 
the contribution of a partition tt to the sum ( 15.1 ip is given by the product 

(5.12) rg, • • -rg, = O L^-W-^WA = 0(n- l ~ 2p+ ^), n -+ oo, 

where |7r e | denotes the number of sets in the partition tt with even cardinality. Since 
\tt°\ < p, the asymptotic contribution of the product (I5.12p is greatest when p = 1, which 
implies tt = U and |vr c | = e(Z). Since the number of terms in the sum (15. lip remains finite 
as n — > oo, by setting p = 1 and combining Eqs. ( I5.12p and ( I5.10P we obtain 

(5.13) n lim = rg(0) = *(/3) 2 (/ - 3 + c(Q) (/ - 4 + e(l)) ^(0) . 
In order to complete the proof of Lemma 15.11 we observe that formula (IB.ip gives 

(5.14) 6{f ) = -L + o (n- 1 ) , n^oo. 
Since when j3 G {1,4}, i(/3) 2 = 4//3, we arrive at 

lim n l+ < 1 ^ (126(f)//?) 1(1 -!)(/- 2)^\ = -^-1(1 -!)(/- 2) 

(5.15) n ^°° J-OP 

x (/_3 + e (/))(/-4 + e (/))^ ) 3 . 
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Combined with the limit (I5.6p . Eq. ( 15.1 5p implies that the asymptotic expansion ( 15. ip is 
valid for every I > 3 by induction. That is, we have 

(5.16) lim n l+e ®- 2 A(l)4l\ = -f3 2 (l + 1)«{2(0). 

n— s>oo 

Finally, inserting these limits into ( 14. 6 p yields the recurrence relation ( 15. 3p . 

5.2. The Solution of the Limiting Recurrence Relation (15. 3p . The limiting recur- 
rence relation ( 15. 3 p is linear and can be solved using the method of generating functions. 
Since either the even or the odd cumulants appear in ( 15. 3p . we introduce 

(5.17a) F W {x):= J2^Mx l 



i=i 



(5.17b) F^(x):=J2 



(2/)! 



;2 (2/-l)r 

Proposition 5.2. Let f3 = 1,4. If at leading order the cumulants of the conductance 
satisfy the recurrence relation (15.31) with initial conditions (15.41) and (15. 5p . the generating 
functions (I5.17P are 

(5.18a) F w {x)= W-a)(a + 6/2 + 2-fl) z/(4/3 2 ) 



2f3 2 1 - x/(4/3 2 



|5 \/ x 



(5.18b) F e W(x) = f | - 1 J f 1 - ^ - " */(4/3 2 

/or -4/3 2 < x < 4(3 2 . 

Proof. In order to simplify the algebra, it is convenient to rescale the cumulants k\ , so 
that the difference equation ( 15. 3 p becomes independent of the parameters a, /3 and 5. 
Define 

(5.19) Si = f$5 K! " )<0) ' ' ^ 3 ' 

where C(a,5, /3) is constant in / and is fixed by the initial conditions. In terms of the 
rescaled cumulants k\ Eq. (I5.3P becomes 

(5.20) - (I + + 1(1 - 1)(4Z - 5)k,_i 

= 3/(/ - - 2)(l - 3 + e(0) (/ - 4 + e(0)K { _ 3 , / > 6. 

This recurrence relation is linear and independent of a, S and f3; therefore, these parameters 
can affect the /cj's only through an overall multiplicative factor, which can be removed by 
an appropriate choice of C(a,8,{3). Then, we introduce the rescaled generating functions 
as 

(5.21a) ^)^^^^ ( ^ ) = f^a B ., 
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The initial conditions (15. 4p and ( 15. 5 p suggest the choice of the constants C Q and C c : 

(5/2 -a) (a + 5/2 + 2- p) 



(5.22a) C (a,6,P) 
(5.22b) C c (a,5,(3) = ( | - 1 



2/3 2 

2 _,v 

Therefore, the rescaled initial conditions are 

(5.23a) k 3 = 2, k 5 = 24, 

3 15 

(5.23b) k a = -, k 6 = —. 

Combining the series expansions in the right-hand sides of (15. 2 lap and (15.21bl) with the 
initial conditions ( I5.23ap and ( 15.23bp . respectively, turns the recurrence relation (I5.20p into 
the remarkably simple ODEs 

dF (v) 

(5.24a) 23/(1 - 3y)(l - V)—^ + (6j/ 2 - 3y + l)F (y) - 3j/(l - 2y) = 0, 

(5.24b) 4(1 - y)*M + 2F e (y) -y = 0, 

dy 

with initial conditions F o (0) = and F e (0) = 0, respectively. These equations can be 
integrated with elementary methods. We have 

y 



(5.25a) F (y) 



i-y 

y 



(5.25b) Fe ( y ) = i_|_yi_ y . 

We give some additional details on this calculation in Appendix [Dl □ 

Remark 5.3. The Taylor series of the right-hand sides of Eqs. ( I5.18a|) and (15.18b|) are 
elementary and prove Theorem 12.11 when I > 3 and k = 0. 

5.3. Joint Cumulants of Conductance and Shot Noise. In order to complete the 
proof of Theorem 12 .![ we need to extend the analysis of the leading order asymptotics of the 
cumulants of G to the joint cumulants of G and P. The strategy follows a similar pattern 
to the approach used for the leading order asymptotics of the conductance cumulants, but 
obviously the technical difficulties increase. 

Lemma 5.4. Let (3 = 1 or /3 = 4 and assume that a is independent of n. Then, for 
(l,k) {(0, 0), (1, 0), (0, 1), (0, 2), (2, 0)} the joint cumulants of the conductance and shot 
noise admit the asymptotic expansion 

oo 

(5.26) k$ ~ n~ v ^ ^(^K 9 , n -> oo, 

g=0 

where u(l,k) = I + k — e(l) and e(l) is defined in (15. 2p . Furthermore, the leading order 
terms nf^ satisfy the double recurrence relation 



(5.27) {21 + 3k + 2)4^(0) - 2(3k\%^) + (^2)^^(0) 

^ 2 (l + k-l-e(l))(l + k-e(l))^U0), 
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with initial conditions given by the sequences 

(5.28a) ftff (0) = {5/2 -a){a + 5/2 + 2- f3) ^ ~ ^ , when I > 3 and odd, 

P\2p) 

and 

(5.28b) «g>(0) = (f " X ) (J 2 ) . when l>±™d even. 

The initial conditions (15.28aj) and f!5.28b|) are obtained by expanding the generating func- 
tion in Proposition 15.21 In the proof of Lemma [5. 14 Sec. 15.11 it was argued that the initial 
conditions k\q are rational functions of n. It immediately follows from Eq. ( I4.16P that 

the joint cumulants are rational functions of n too and therefore have an asymptotic 
expansion in powers of n~ l as n — > oo. 

The precise form of the expansion (15.26 p . as well as the value of the exponent v(l,k), 
will be determined by mathematical induction on k. First one must deal with a number 
of base cases. Note that if k = 0, Eq. (I5.26P reduces to our earlier result (15. 1ft : if k — 1, 2, 
Eq. (I5.26P follows easily from the recurrence (I4.16P provided 1^0. If k = 3, we again 
use (14.161) to obtain (I5.26p . this time for any I > 0. In the following we shall proceed by 
induction on k to prove (I5.26P for any I > and k > 3. We discuss the left- and right-hand 
sides of Eq. (I5.27P separately. 

5.3.1. Proof of the left-hand side of Eq. ( I5.27P . By rescaling and taking the limit of the 
left-hand side of Eq. (I4.16P we obtain 

(5.29) lim n l+k + 1 -^ (%KC fe -x " " 2(a + 5/2 + nfi + 2 - P)k\% ± 

+(2/ + 3A + 2) K ;« +1 + 6,<«* : g(*: 1 )t(')' 

i=0 v / i=0 VJ/ 

= (2/ + 3k + 2)^(0) - 2/3^ 2ifc (0) + {k/2)K\% k _M- 

In order to prove this limit, one first has to verify that it holds for the initial condi- 
tions (I5.28P by direct substitution. Then, by assuming the expansion (15.261) . we see that 
the only leading order contributions from the double sum in (15.291) come from the choice 
of indices = (0,0) and (J,i) = (l,k — 1). Elementary manipulations involving the 
limiting value of the variance of the conductance give the right-hand side of (I5.29p . 

5.3.2. Proof of the right-hand side of Eq. (15.271) . As in the proof of Lemma ISTTl deriving the 
right-hand side of the recurrence relation (I5.27P presents more technical difficulties than 
the left-hand side. The reason is that the finite-n recurrence (14.161) involves the coefficients 
lJf k defined in terms of the auxiliary recurrence (14.171) involving the two dimensional 
sequence of reduced cumulants 

(5.30) rg? := 



1 h 'j+2,i' % l-j+2,k-i-l 



l.k 



2k\ ) . 



Since we assume that the parameters a and 5 are independent of n, we can insert the 
expansion (I5.26P into (15.301) to obtain 



oo 

(5.31) rg^n-^-^rJfC?^, n oo. 

g=0 
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The coefficients in this series may be computed in terms of those in the expansion (j5.26p . 
At leading order in n, we find 

= z(f3f{l + k-e(l)){l + k + l-e(l))^(0). 
Let us write the coefficient (if^ as a sum over partitions, 

( 5 - 32 ) rffc = XI II r |L(B l )|,|Af(B l )|- 

7T BieiT 

The sum varies over all partitions n = {Bi, . . . , B p }, 1 < p < l+k, of a set U = {1, . . . , l+k} 
into disjoint subsets Bi. For each partition ir, and for each subset Bi, we further divide Bi 
into the two subsets L(Bi) := {j G Bj : j < 1} and M(Bi) := Bi \ L(Bi). We argue that 
there is only one partition tc* which gives the desired asymptotic contribution as n — > oo. 

The contribution of a partition n = {B±, . . . , B p } to the sum (I5.32p is given by the 
product 

r \L(B 1 )\,\M(B 1 )\ X X r \L(B p )\,\M(B p )\ ~ U { n 

(5.33) = O (n-^=^-< L ^ +2 A 

= 0(rT l - k - 2p+ \<Y 

where denotes the number of subsets Bi in the partition 7r such that L(Bi) has even 
cardinality As with the conductance, the dominant contribution to (15.331) corresponds to 
p — 1, i.e. when the partition tt = n* = U. In this case, we have L(U) = {1, . . . , 1} and 
M(U) — {I + 1, . . . , I + k}, and so |7rfJ = e(Z). Therefore, we find that 

(5.34) lim n l+k+1 -<%^) = z(/3)2(/ + k _ 1 _ e(/) ) (j + k _ e (l)) ^(0) . 
Using the leading order of (15.141) and the equality i(/3) 2 = 4//3 we arrive at 

(5.35) lim n'^-O (^k^-\ = JL(l + k - 1 - e(Z)) (I + k - £(0)^(0). 



Combining Eq. ( I5.35j) with the limit (I5.29P implies that the asymptotic expansion ( 153 
is valid for the stated values of / and k. Inserting these limits into the finite-n recurrence 
(I4.16P finally completes the proof of Lemma 15.41 

5.4. The Solution of the Limiting Recurrence Relation. We now complete the proof 
of Theorem 12. II by completing the derivation of Eqs. (I2.4ap and (I2.4bp . 

Firstly, we remove the dependence on a, /3 and S from the recurrence relation (I5.27P and 
its initial conditions (I5.28P by rescaling the joint cumulants k\^. We write 

f5 36) k - ^IM^1 K W( ) 

(5-36) K, k - c{atM «,,*(0). 

As discussed in Sec. I5.2[ since Eq. (I5.27P is linear, it is not affected by the constants 
C(a, (3, 5), which are fixed by the initial conditions (I5.28a|) and ( 15.28b|) ; therefore, for con- 



venience, we choose the same as in Eqs. (I5.22a|) and f!5.22bl) . Substituting (I5.36P into (I5.27P 
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yields 

(5.37) (2/ + 3k + 2)R, l>k+1 + 2kk l+2 , k -i - 2k i+2ik 

= 3k(l + k - 1 - e(Z)) (l + k — e(l))k l>k ^. 
We treat separately the cases / even and / odd. Replacing Z with 21 gives 

(5.38) (4Z + 3fc + 2)re wl + 2kK 2 i+2,k-i ~ 2k 2l+2 , k = 3k(2l + k- 2) (21 + k - l)R 2 i,k-i- 
We need to solve this difference equation subject to the initial conditions 

(5.39) ^W-W 



A 1 V / 

We claim that the unique solution is given by 



(5 , 0) ^ = ( _ 1)4 PLt*^!g(^f)g 

This formula is proved in Appendix O using the properties of Jacobi polynomials. 
For the odd cumulants, replacing / with 21 — 1 in (I5.37P leads to 

(5.41) (4/ + 3k)R 2l ^ k+1 + 2fc«2i+i,fc-i - 2k 2l+x>k = 3k(2l + k- 2) (21 + k - l)R 2 i-i, k -i- 
We need to solve this recurrence subject to the initial condition 

(5.42) Ka_ 1>0 = (2Z-2)! 

It is straightforward to prove by direct substitution that the unique solution is 

(5.43) K2i-i lfc = (k + 21-2)1. 



This complete the proof of Theorem 12.11 



5.5. Higher Order Asymptotics. As we discussed in Remark I2.2[ the cumulants of the 
shot noise ( I2.6a|) of odd index appear at sub-leading order in the large-n expansion (|5.26p . 
Therefore, in order to prove fl2.6a[) . we need to calculate KqI(1) for k odd. The purpose of 
this section is to give a flavour of the difficulties one encounters when trying to compute 
next to leading order terms of the cumulants using our approach. 

In order to simplify the analysis, we consider the case with j3 = 1, a = — 1/2 and, 5 = 0, 
i.e. the scattering matrix belongs to the COE and both leads carry the same number of 
channels. 

Inserting the asymptotic expansion (I5.26P into the recurrence relation (I4.16p . setting 
I = and equating terms of order n~ k for k odd yields 

(*/2)«gL(l) + (3* + 2)«£Ul) - 2 «S(1) " (3/16)fc 2 (A; - 1)k£L(1) 

fc-2 



+ 6k T,( k i ^S-i-iW^W- (3^/4)^(0) -«g(0). 



From this equation we can see that the cumulants k^\{V) are expressed in terms of the 

mixed cumulants K^fc-iC-O anc ^ K 2kQ) m Thus, it becomes clear that to compute the KQ k (l)'8 

one requires the knowledge of d^ik^X) ^ or eacn ^ 

Define gi jk '■= K 2Z Then, at next to leading order Eq. (I4.16P becomes 
(k/2)g l+ltk _ x - 2g l+1>k + (4Z + 3k + 2)g ltk+1 - (3A;/16)(Z + k — 1)(Z + fc)0i,k-i = 6(Z, k), 

where 8(Z, k) is a complicated function of k 2 ] \(0). 
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Using MAPLE we were able to surmise that for I > 1 the quantities 

. 4^0) 4«^ 2ifc (i) 

l ' k (21 + k-iy. (2i + k + iy. 

can be extracted from the generating function 

h Uk = Res (Res (x~ l - 4 y- k+1 F(x, y/x)y 2 )) , 

where 

/ rv V¥P(x,y') j A 

- 8y + A)y) m 

The function P(x, y) is the polynomial 

P{x, y) = -512 + 256x - 20x 2 + 64y V - 128yx + y V - 16y 2 x 3 + 3yx 3 - 8yx 2 

We remark that the integral in Eq. ( I5.44j) can be computed in closed form, but it is 
rather long and cumbersome to be reported here. 

6. The Wigner Delay Time 

The analysis of the random variable ( II. 6p follows the same strategy that we adopted 
to study the conductance. The main difference is that the exponent (11.111) depends on n 
explicitly, which means that the analysis of the asymptotic behaviour of the cumulants is 
more complicated. A direct study of the integral (12.151) would require the solution of a 
delicate double scaling problem. In addition, the linear homogeneous recursion relations 
obtained in the previous section must be replaced with non-linear recursions that are harder 
to solve. Another difficulty is the singular nature of the weighting function (13.31) . which 
implies that the cumulant generating function is not analytic and only a finite number of 
moments exist. 



6.1. The Differential-Difference Equation. The proof of Theorem 12.91 is entirely anal- 
ogous to that of Theorem 12 .51 Therefore, we shall only outline the relevant steps. 
Recall that the r-function associated to the Wigner delay time is 



— / 77 Pz(Xj) ex P ( ^ tj\) j 71 l^fc _ XjfdXi ■ ■ -dX n , 

n - JR™ =1 V i=1 / Ki<k<n 



(6.1) r^(t;z 

mi,-,, 

j=l \i=l / l<j<k<n 

where p z (X) was defined in Eq. (13.31) . It is related to the moment generating function by 
the equation 

MW(^ y /2) = -^r(t;/3y/2) . 

tn n t=o 

Because of its universal usage and since there is no risk of ambiguity, we adopt the same 
notation to indicate the r-function that we used for the conductance and for the shot noise. 

This r-function satisfies the KP equation (13. 9p or Pfaff-KP equation (13. 8p depending on 
the value of (3. The definition (16. ip gives 

(6.2) ^W {t] z) = ^(t- i z). 
Thus, the derivative 

Kn [) ~ dz 
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can be expressed in terms of the projection at t = of the right-hand side of Eq. (16. 2p . It 
remains to express the partial derivatives 







2 



d 



,2 



0l - 2 ' - . „ > , and -r^r log (t; z) 



t=0 dt x dt z 



t=0 



(6.3) ^logr^(t;z 

in terms of the derivative of £n (2). The extra relations needed for this purpose are 
provided by the Virasoro constraints. 

The only Virasoro operators needed for these calculations are V-i and Vo- They are 
expressed in terms of the parameters (a»)^o anc ^ (^)£o defined in Eq. (I3.10p . which are 
derived from the relation 

d\ogp z {\) _ z\ 2 - bX + 1 
dX ~ X 2 ' 

This equation gives 

a 2 = 1, b = — 1, bi = b, b 2 = -z, 

and otherwise = and ft, = 0. Finally, a lengthy but straightforward calculation 
analogous to that in Appendix |A] leads to the ODE (12.2 



6.2. Finite n Recursion Relation. Let us introduce the reduced cumulants 

-V) 



(6.4) tf> := ^ + K* 

n—tn—iyp) 



-2Kf\ 

n— 5>n+i(/3) 



where the K\ s are the coefficients of the expansion (12.211) . as well as the reduced moments 

3=0 V J J 



We also define 

(6.5a) Ail) := r][^l(l - 1)(/ - 2) + 01(31 - 2) - w(w - 4 + 2(3)1 - w(w + 2 - /3), 

(6.5b) := (j - i) (6 V [% -i-l)i + + 4/3z) , 

where a; := 6 — 2 — (3n + (3. Recall also that q := |_& — 2 — /3(n — 1)J . 

Lemma 6.1. The cumulants of the Wigner delay time are the unique solution of the 
recurrence relation 



11 'I 



(6.6) a{i)k^\ +J2() B ^ o*3*ffi + - m 



(3) 



i=0 

= (12df / 0)1(1 -l)(l- 2)09/2) VS, 3 < / < 9 - 1, 
with initial conditions given by Eqs. (12.321) . 

The proof is entirely analogous to that of Lemma 14.11 Inserting the expansion (I2.2ip 
into (12.201) and equating the first three powers of z gives a system of linear equations whose 
solutions are the initial conditions (12.321) . Equating higher powers lead to (I6.6p . 

Remark 6.2. Equation (16.41) depends implicitly on the exponent b that appears in p z (X). 
Since the Pfaff-KP equation (13. 8p was derived with the assumption that the weighting 
function in the integral (13. ip should be the same for any dimension n of the r-function, 
the reduced cumulants are well defined only if all the Ki's in the right-hand side are 
computed for fixed b. Once a solution (ifj )? =1 of the difference equation (16.61) is found, 
the exponent b can be assigned the value (II. lip to recover the cumulants. 
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The implicit dependence of the reduced cumulants (16.41) on the exponent (11. lip means 
that the asymptotic analysis of the cumulants of 7\y f° r — 1>4 involves a double scaling 
limit; so does a direct study of the asymptotic behaviour of integral representation (I2.15P 
of M n (z). These problems go beyond the scope of this paper and remain open. In the next 
section we shall study the asymptotic limit of the cumulants with = 2 symmetry classes. 
In this case the right-hand side of the recurrence relation (16.61) is zero and a direct study 
of the limit n — > oo presents less challenges. 

6.3. The limiting Recurrence Relation: /3 = 2. The proof of Theorem 12.121 combines 
the approaches in Lemma 15.11 and Proposition 15. 2} which provide the weak localization 
corrections for the cumulants of the conductance. The main difference is that the limiting 
recurrence relation is not linear. However, it can still be solved using the method of 
generating functions. 

Equations (16.51) and (16. 6p implies that the cumulants are rational function of n. A direct 
calculation shows that admits the asymptotic expansion 

oo 

(6.7) ~ n 2 - 21 K^\p)n- p , n -> oo, I > 0. 

p=0 

This formula is correct for all the symmetry classes G {1,2,4}. Note that the exponent 
in front of the sum does not contain a term analogous to z/(Z); therefore, we do not see the 
"staircase" asymptotic behaviour which characterized the cumulants of conductance and 
shot noise (see Remark 12.21) . 

When = 2, our fmite-n recurrence relation (16.61) simplifies considerably: 

(6.8) (/ + l)(n 2 - l 2 )Ki 2 \ = 21(21 - 1)K {2) + 2 f> + 1) ( l \ (I - ifK^K^, 

i=0 ^ ' 

where we made use of the identity 

(6.9) ( ■) ( l ~ *)( 6 ( Z - * - 1)* + 2 + tyxi+txt-i = 2 ^(3z + 1) Q (I - ifx i+1 x^i, 

i=0 ^ ' i=0 ^ ' 

which follows from equating coefficients of the indeterminate Inserting the as- 

ymptotic expansion (16.71) into (I6.8P and taking the leading order terms gives the limiting 
recurrence (I2.28p . after identifying the new variable 

K, (2) (0) 

(6.10) p, = (iAi. 



Difference equations similar to (I2.28P have been studied in [42l . |48| using formal power 
series solutions of Painleve I. In a similar spirit Garoufalidis et al. [39j used the Riemann- 
Hilbert approach to study asymptotics of non-linear recurrence relations. 

Recall that 



F(z)=J2Pi 



i=i 



This generating function F(z) turns the difference equation (12.281) into the non-linear ODE 

(6.11) 2F{z) + F'(z) - AzF' - QzF'(z) 2 + 4F(z)F'(z) = 1, 

with initial condition -F(O) = 0. The constant in the right-hand side is fixed by the 
condition that the first cumulant is normalized to one. Equation (I6.1ip is of d'Alembert 
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type and can be integrated. We have 

(6.12) F(z) = 3/2 - 3/2Q(z) + 2zQ(z) 3 + 3zQ(z) 2 + 2zQ(z), 
where Q(z) is the unique solution of the quartic equation. 

(6.13) 4zft 4 + 8zQ 3 + (4z - 3)Q 2 + 2Q + 1 = 0, 

with tt(0) = 1. 

Solving (I6.13P for z in terms of Q gives 

Since, as a function of fl, z(Q) is analytic at 1, z{l) = and z'{l) ^ 0, we can apply the 
Lagrange inversion formula and obtain the formal power series 



(6 ,5) F= 1+ 5>* 



fc=l \ \ / / fc=1 



where 



(6.16) 




fe-l-i 



The right-hand side follows from computing the derivatives and rewriting the resulting 
factorials in terms of binomial coefficients. This completes the proof of Theorem 12.121 

Since the factor ( 2 ^) is divisible by k, we see that the Cfc's are integers. It follows that 
the Taylor coefficients of the generating function f)6.12p are integers too, and therefore the 
solution to fT2T28]) . 

A natural question would be to ask whether the cumulants of the Wigner delay time for 
the symmetry classes (3 = 1 and (5 = 4 are integer too. In Tabled we list the leading order 
coefficients of the first eight cumulants, which we computed numerically, for (3 e {1,2,4}. 
We immediately see that when (3 = 4 this is not the case is not the case; when ft = 1 it 
remains an open question. 



Table 2. The limiting cumulants of the Wigner delay time for each ft G {1,2,4} 



1 


^ (1) (o) 


K^(0) 


^ 4) (o) 


1 


1 


1 


1 


2 


4 


2 


1 


3 


96 


4 


6 


4 


5088 


636 


159/2 


5 


437760 


27360 


1710 


6 


53038080 


1657440 


51795 


7 


8353013760 


130515840 


2039310 


8 


1625430159360 


12698673120 


396833535/4 
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7. Conclusions and Open Problems 

In this paper we study several properties of the cumulants and of their generating func- 
tions of the conductance, shot noise and Wigner delay time in ballistic quantum dots for 
all the symmetry classes /3 G {1,2,4}. We focus on the behaviour for finite numbers of 
quantum channels n and at leading order as n — > oo. The number of channels can be 
interpreted as the semiclassical parameter. 

The main results are the limiting formulae as n — > oo of the cumulants given in The- 
orems |2J] and |2T2j as well as the differential equations (12. lip . (12.13p and (12. 20 p . which 
are satisfied by the cumulant generating functions at finite n. Theorem 12.11 provides the 
weak localization corrections of the mixed cumulants of the conductance and shot noise 
for the symmetry classes ft = 1,4. Theorem 12.121 concerns the leading order asymptotics 
of the cumulants of the Wigner delay time when (3 = 2. The differential equations (12. lip . 
(I2.13P and (I2.20p apply to all symmetry classes (3 6 {1,2,4}. When (3 = 2 we show that 
the cumulant generating function of the Wigner delay time can be expressed in terms of 
the Painleve III' transcendent and discover an interesting identity between the probabil- 
ity distributions of the conductance and shot noise. Finally, we derive finite n as well as 
asymptotic recurrence relations that are very effective to compute cumulants of arbitrary 
order. 

The results for the conductance and shot noise apply not only to the Wigner-Dyson 
symmetry classes, but also to those introduced by Zirnbauer [§8j and Altland and Zirn- 
bauer [§, [Io[ , which characterize superconducting quantum dots (also known as Andreev 
quantum dots). However, the formulae for the Wigner delay time are limited to the Wigner- 
Dyson symmetries. 

A number of conjectures proposed by Khoruzhenko et al. [49[ are proved as particular 
cases of Theorem 12.11 They concern the weak localization corrections of all the cumulants 
of the conductance (see Eqs. (12.5aj) and ( 12.5b|) ) as well those of the even cumulants of the 
shot noise (see Eq. (12.6b|) ). In the same article they conjectured the leading order ( |2.6aj) 
of the odd cumulants of the shot noise. Unfortunately, the latter are of subleading order 
to our formulae; their proof remains an open problem. 

Our study is far from being exhaustive and several interesting questions remain unsolved. 
Our approach is based on the integrable theory of r-functions that are deformations of 
certain matrix integrals. These techniques were mainly developed by Adler and van Mo- 
erbeke (H^4s|] and Adler et al. ^M^. While these ideas appear to be very effective to study 
the leading order asymptotics of the cumulants, going beyond the leading order seems a 
very difficult task. We discuss this aspect in Sec. 15.51 

Another challenging project would be to compute the leading order asymptotics of the 
cumulants of the Wigner delay time when (3 = 1,4. Such computations would involve 
a double-scaling limit that appears to be intractable using the techniques in this paper. 
Interestingly, we discover that at leading order the cumulants of the Wigner delay time are 
integers when (3 = 2, but are not when /3 = 4; when (3 = 1 it remains an open problem. 

There is some evidence that the limiting formulae in Theorem 12.11 hold for any real 
(3 > (see Remark 14. 5p . i.e. for the /3-ensembles, which are of significant interest in 
many applications of modern RMT. Generalizing the techniques applied in this paper to 
arbitrary (3 > poses some fundamental problems, namely it is not known whether there 
exist analogues of Eqs. ( 13. 9 p and (13. 5p for other values of (3. 

We conclude by pointing out that Dahlhaus et al. [28[ have suggested that the Andreev 
quantum dots of symmetry types Dili and CI - in the Cartan classification - - could 
be realized experimentally using graphene. This opens the possibility that some of our 
formulae could be verified experimentally. 
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Appendix A. Proof of the Differential Equation (I2.13p 



The derivations of the ODEs in Theorems I2.4[ 12.51 and 12.91 are based on the same ap- 
proach. The projections at the origin of some the partial derivatives of the r-function in 
the KP and Pfaff-KP equations (I3.9P and (13. 8p are expressed in terms of the derivatives of 
the moments generating function using the Virasoro constraints. The calculations are in- 
volved, but standard in the theory of integrable systems. Here we discuss the derivation of 
the ODE (I2.13p . which is satisfied by the cumulant generating function of the conductance. 
In this case we need to write the projections (}4.4j) in terms of and its derivatives. 

For ease of notation, let us define 

0l : = (n/3 + 2-/3), o 2 := (n/3 + 3 - 3/3/2), c := an + n({n - l)p/2 + l). 
Recall that 

a^(z)=logMf(z) = logr n (t;z)\ t=0 

and 



<9 fc logr„(t; z) 



dt\ 



-1) 



t=o 



d k a^\z) 
dz k 



Setting w = in (13. 16a|) and (I3.16bp . we arrive at the following Virasoro constraints: 



(A.la) 



V_i logr n (t; z) 



,i=i 







d 



dt i+1 dti 



+{a + Z + S/2 + oi) 



_d_ 



_d_ 



c ) logr n (t; z) = 



and 
(A.lb) 



V logr n (t;z) 
^(t) 



i=l 











dt 



i+2 



dt 



i+l 

d 



+ 



2dt\ 
d 



Wi + a. 



_d_ 



+ (a + z + 5/2 + o 2 )— - z— ) logr n (t; z) 



+ 



(3 ( d\ogT n (t;z) 



0. 



2 V dt! 

Setting t = in (lA.lap and (lA.lbl) . respectively, leads to the relations 



91ogr n (t; z) 



dt 2 

91ogr n (t; z) 



1 



-(c+(a + z + 8/2 + 0l )af)., 



t=o 



t=o 



z \2 n 2 { n ] 
+ (oi + a)af + (a + z + 6/2 + o 2 )X x ) . 
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Differentiating both (I A. lap and (lA.lbj) with respect to ti and setting t = yields two 
further relations 



X 4 :~- 



d 2 logr n (t; z) 



dhdt 2 
d 2 \ogr n (t;z) 



dtidtz 



t=0 



t=0 



l(ap + X 1 + (a + z + d+o 1 )ap'), 
- (01 + a)aW" + (a + z + 5/2 + o 2 )) . 



Taking the partial derivatives of the constraint ( lA.lal) with respect to t 2 and setting t = 
gives 



d 2 logr„(t:z) 



dt\ 



■ (2X 2 - 2X 1 + (a + 5/2 + z + 0l )X 3 ) . 



t=o 



Therefore, each Xi, i — 1, ... ,5, can be expressed solely in terms of known parameters and 
derivatives of Finally, evaluating the Pfaff-KP and KP equations (I3.8P and (13. 9p 

at t = and inserting X 4 and X 5 yields (I2.13p . 



Appendix B. The Constants 6n /3) and dtP 

B.l. Conductance and Shot Noise. The constant bn \ as introduced in Eq. (I2.9p . is 

defined in terms of Selberg's integral, whose explicit expression is given in Eq. (jl.Sp . 

The ratio of Selberg's integrals in (I2.9p involves a large number of cancellations of Gamma 
functions. Straightforwards algebra gives 



(B.la) &W :- 



I 1\ (!) (!) 

n(n - 1) c^ 2 <j 2 



(n+l)(n + 2) ( C W) 2 

n(n - l)(2a + n)(2a + n + 1)(5 + n)(5 + n + 1) 
16(5/2 + a + n)(S/2 + a + n + l) 2 (5/2 + a + n + 2) 
(5 + 2a + n + 1)(5 + 2a + n + 2) 



x 



while for — 4 
(B.lb) 6^ := 



(5 + 2a + 2n - 1)(5 + 2a + 2n + 1) 2 (5 + 2a + 2n + 3) : 



(4) (4) 
71 C n-l C n+l 



" + 1 (ci 4) ) 2 

2n(2n + l)(a + 2n)(a + 2n - l)(5/2 + 2n)(5/2 + 2n - 1) 
(5/2 + a + 4n)(5/2 + a + An - 2) 2 (5/2 + a + An - 4) 
(5/2 + a + 2n - l)(5/2 + a + 2n - 2) 



x 



(5/2 + a + 4n + l)(5/2 + a + An - l) 2 (5/2 + a + An - 3) ' 



B.2. Wigner Delay Time. The parameter dn , see Eq. ( I2.18p . has the same structure as 
b^n \ with the difference that Selberg's integral is replaced by the normalization constant 
, which appears in the definition of the partition function (I2.15p . Similar calculations 
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to those leading to Eqs. (IB. lap and ( IB. lb|) give 



n(n — V 



(i) (i) 

m n-2 m n+2 



(n + l)(n + 2) ( m «) 2 

n(n - 1)(26 - 2 - n)(26 - 1 - n) 
(b - n)(2b -2n + l)(b-n- 2) (26 - 2n - 3) (6 - 1 - n) 2 (2b - 2n - l) s 



and 



(4) (4) 
n m n-l m n+l 



71+1 



(4)\2 



K 4j ) 



2n(2n + 1)(6 + 2 - 2n)(6 + 1 - 2n) 



(b + 3 - 4n)(b + 1 - 4n) 2 (6 + 2 - 4n) 2 (6 - 1 - 4n)(6 - 4n)(6 - 4n + 4)' 

B.3. Dualities. One easily sees that the constants &n and display the following 
dualities between the symmetry classes /3 = 4 and — 1: 



(B.2a) 
(B.2b) 



6^ = &« 

"n n 

16dl 4 ) = dW 



<5— ► — 5 / 2 , a — > — a / 2 , n — > — 2 n 



6-+- 6/2,n->— 2n 

The duality^ (IB.2al) . and other similar relations, were also mentioned by Adler and van 
Moerbeke 



Appendix C. Proof of Eq. (I5.40p 

We now prove that (15.401) is the unique solution of the initial value problem (15.381) with 
initial conditions (15.391) . 

First observe that the right-hand side of Eq. (15.401) is expressible in terms of Jacobi 
polynomials using the identity 

(CD t( 2 ^f)Q(-2)--m^-(-3). 



Recall that the Jacobi Polynomials can be defined by 

f rt + a\ (n + b\ f x — 1\ n ~ 3 ( x + 1 x 1 

J + K 



3=0 



J 



The identity (IC.ip is a consequence of this formula. 

Inserting the right-hand side of (IC.ip into (15.401) and substituting K2i,k into (I5.38P gives 



(C.2) (4/ + 3k + 2){k + l)P^ k ' 3/2 \-3) + {21 + !)(* + I + l)P fc ^' 



(Z+l,-fc+l/2) . 



-3) 



+ (21 + k){2l + l)P, 



(l+l-k-l/2) 



(_3)_3(fc + Z )(fc + J + l)p£- 



-3) = 0. 



This identity can be proved using well known properties of the Jacobi polynomials. Firstly, 
we identify the parameters of each Jacobi polynomial in (1C.2I) using the identities (see, 
e.g., [1], Eqs. 22.7.18 and 22.7.19) 



{2n + a + b)P { n a ^ 1 \x) 
{2n + a + b)Pi a ^' b \x) 



[n 



71 



+ a + b)P^ b \x) + {n + a)P^< 



x), 



+ a + b)P^ b \x) - (n + b)P^{x 
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Then, we repeatedly apply the three-term recurrence relation (see, e.g., Eq. 22.7.1) 

p( a,b) ( \ = 2{n + l){a + b + n + l){2n + a + b)P^(x) 
%n {) {a + b + 2n)(a + b + 2n + l){a + b + 2n + 2) 

{2n + a + b + l){b 2 - a 2 )Pi a ' b \x) + 2{a + n)(b + n){a + b + 2n + 2)P%\{x) 
+ {a + b + 2n){a + b + 2n + l){a + b + 2n + 2) 



We used a similar technique to prove equalities such as ( \C.2\i in a previous publication [63 
(see, e.g., Proposition 3.10). 

Appendix D. The Solution of the Recurrence Relation (15. 3p 

The recurrence relations (15.31) and (12.281) can be solved using the method of generating 
functions. We now give a flavour of the main ideas involved by discussing in some detail 
the example of Eq. (15.31) . 

The change of variables (15.191) turns Eq. (15.31) into (15.201) . Let us first consider the even 
cumulants. By replacing / with 21 — 1 in (15.201) . we obtain 

(D. 1) - 2lk 2l + {21 - 1) {21 - 2) (8/ - 9)k 2l _ 2 = 3(2/ - 1) {21 - 2) {21 - 3) {21 - 4) {21 - 5)k 2 i-4, 

which we need to solve subject to the initial conditions k 4 = 3/4 and k e = 15/2. Introduc- 
ing the sequence pi = k 2L /{2l — 1)!, / = 2, 3, . . . , the recurrence (ID.ip becomes 

(D.2) - 2l Vl + {81 - 9)p z _! = 3(2/ - 5)p,_ 2j 

with initial conditions p 2 = 1/8 and ps = 1/16. 

The difference equation (1D.2j) is transformed into a differential equation for the gener- 



ating function F e {y), defined in (15.21bl) . by multiplying both side by y l and summing over 
/ from 4 to oo. For example, we have 

oo , oo , 

i=a y 1=4 y 

Repeating this procedure for each term in (1D.2|) and simplifying leads to 

dy 



dF 

(D.3) y{3y - 1) I y - 2F e + A(y - 1)— ^ I = 



with initial condition -F e (0) = 

For any y in the punctured disk {y G C | < \y\ < 1/3}, we divide both sides of (ID.3j) 



by y{3y — 1) and solve the resulting first order linear equation by the method of integrating 
factors. The solution is 

y 



F e {y) = l-|-Vl-y. 

The odd cumulants can be found in a similar fashion. The variables pi = k 2 i + i/{2l)\ 
satisfy the recurrence 

(D.4) - {21 + l) Pl + (8/ - 5)p,_i = 6(Z - l)pi- 2 , 

with initial conditions p\ — 1 and p 2 = 1. Then, one can show that the generating function 
F Q {y), defined in Eq. (I5.21a|) satisfies the linear ODE 

2y(l - 3y)(l - y)^^- + {6y 2 -Sy + l)F {y) - 3y(l - 2y) = 0, 
dy 
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with initial condition -F o (0) = 0. This equation can be solved by the method of integrating 
factors. The solution is 

F (y) = ~, ~~ • 
i - y 

Alternatively, one easily verifies by direct substitution that pi — 1 for all / > solve 
Eq. flDU). 
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